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ABSTRACT 


Miles'  theory  is  unsatisfactory  for  the  prediction  of  ocean  wave 
growth  under  wind  action.  Hasselmann's  nonlinear  wave-wave  inter- 
actions theory  may  offer  another  physical  mechanism  for  constructing 
a reasonable  wave  prediction  model.  An  experimental  program  was 
designed  to  investigate  this  possibility  in  the  context  of  the  momem- 
tum  and  energy  transfer  from  wind  to  waves,  nonlinear  wave-wave 
interactions  and  white-capping  dissipation.  The  wind  field  (mean 
velocity  profiles  and  fluctuating  pressure  and  velocity  components) 
and  wave  field  (wave  height)  were  measured  simultaneously  in  a fixed 
reference  frame  and  as  a function  of  fetch  along  the  Stanford  wind- 
wave  channel  under  the  conditions  of  steady  wind  and  stationary  wave 
spectrum.  All  the  data  were  obtained  5 mm  above  the  highest  point  of 
the  wind  waves  for  five  stations  (3  m apart  on  average)  and  at  three 
wind  speeds  (7.09,  8.01  and  8.88  m/sec).  The  wave  height,  fluctuating 
pressure  and  velocity  components  were  measured  by  a capacitance  wave 
height  gauge,  a crystal  pressure  transducer  and  a cross  hot  film  probe, 
respectively. 

The  normalized  parameters  of  a wave  field  such  as  spectral  peak 

A A 

frequency  f , total  energy  e and  the  Phillips  equilibrium  constant  a 
m 

of  the  wave  spectrum  were  found  to  be  fetch  and  wind  speed  dependent. 
The  power  spectra  of  turbulent  velocity  components  in  the  inertial  sub- 
range were  consistent  with  Kolmogorov  -5/3  power  law.  The  small  scale 
structure  of  the  power  spectra  of  turbulent  velocities  seemed  to  be 
wind  speed  dependent  but  fetch  independent. 
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The  energy  and  momentum  transfer  from  wind  to  waves  were  confirmed 
to  be  dominated  by  the  normal  pressure  acting  on  the  wave  surface.  In 
addition,  the  oscillating  wave-associated  turbulent  stresses  accounted 
for  about  30%  of  the  total  momentum  transferred  into  waves,  but  made  a 
negligible  contribution  to  the  energy  transfer  from  wind  to  waves. 

The  average  ratio  of  momentum  transferred  into  waves  over  the  total  mo- 
mentum transferred  across  the  wind-wave  interface  was  found  to  be  about 
0.65.  The  remaining  35%  of  the  total  momentum  transfer  goes  directly 
to  currents. 

Based  on  the  experimental  results,  Hasselmann’s  nonlinear  wave-wave 
interaction  theory  appears  to  be  valid.  Barnett's  approximate  parame- 
tric equation  for  calculating  the  energy  transfer  of  nonlinear  wave- 
wave  interaction  and  Hasselmann's  white-capping  dissipation  model  were 
also  verified  and  appeared  to  be  applicable  in  the  relatively  low  and 
intermediate  frequency  region  of  a wave  spectrum  for  a normalized  fetch 
range  of  100  < x < 500. 

A method  of  modifying  Barnett's  parametric  equation  representing 
the  nonlinear  wave-wave  interaction  in  the  relatively  high  frequency 
region  was  suggested.  Based  on  the  results  of  overall  energy  balance 
of  a gravity  wind  wave  spectrum,  the  nonlinear  wave-wave  interaction 
mechanism  was  confirmed  to  play  a dominant  role  in  the  energy  transfer 
processes  after  the  wave  spectrum  is  generated. 
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Fourier  transform  of  x(t)  and  y(t) 


Greek  symbols 


a 


Phillips  equilibrium  constant  of  wave  spectrum 
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Phillips  equilibrium  constant  of  wave  spectrum 
(=  lb’ir^'o) 

angle  between  k2  and  (k^  - k^)  (radian) 

ratio  of  momentum  transferred  into  waves  over  total 
momentum  transfer  across  the  wind-wave  interface 
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ratio  of  energy  transferred  into  waves  over  total 
energy  transfer  across  the  wind-wave  interface 

damping  factor  in  Hasselmann’s  white-capping  dissipation 
model  (-  d (2irf)2) 

one  and  two  dimensional  functions  in  Barnett’s  parame- 
tric equation  defined  by  equations  (2.18)  and  (2.15) 

pressure  difference  from  the  pitot-static  tube  measured 
across  the  Pace  pressure  transducer  (psi) 

total  wave  energy  (»  jE(f)  df)  (cm2) 

normalized  total  wave  energy  (=  eg2/u^Q) 

normalized  total  wave  energy  (=  eg2/u£) 

water  surface  (wave)  elevation  (cm) 

Fourier  transform  of  n 


V 


xv  i 


conjugate  Fourier  component  of  TT(f) 

directional  angle  measured  from  north  (radian) 

mean  direction  parameter  of  Barnett's  parametric 
equation  (radian) 

phase  angle  between  variables  x and  y (degree) 

von  Karman  constant  (=  0.4) 

a linear  factor  between  X(f)  and  fl(f) 

kinematic  viscosity  (cm2/sec) 

density  of  air  (gm/cm3) 

density  of  water  (gm/cm3) 

one  and  two  dimensional  functions  in  Barnett's  para- 
metric equation  defined  by  equations  (2.19)  and  (2.16) 

mean  momentum  transfer  due  to  white-capping  dissipa- 
tion, air- input  and  nonlinear  wave-wave  interaction  in 
the  high  frequency  region  (dyne/cm2) 

mean  momentum  transfer  from  wind  to  waves  due  to  normal 
pressure  acting  on  wave  surface  (dyne/cm2) 

phase  angle  of  the  transfer  function  H(f)  of  the 
modified  system  of  crystal  pressure  transducer  (degree) 

radiant  frequency  (»  2irf)  (radian  - Hz) 


I . INTRODUCTION 


1.1  History  and  Motivation 

In  recent  years  many  scientists  from  varied  backgrounds  have  been 
working  cn  problems  of  wind  wave  generation.  Although  important  ad- 
vances have  been  made  in  the  wind  wave  field  during  the  last  twenty 
years,  a basic  understanding  of  the  dynamical  processes  occurring  at 
the  air-sea  interface  has  not  been  achieved.  The  difficulties  which 
impede  the  progress  of  the  air-sea  interaction  problem  are  the  lack  of 
fundamental  knowledge  of  the  nature  of  interaction  between  the  turbu- 
lence in  the  air  stream  and  the  perturbed  water  surface,  and  the 
complex  and  elusive  nature  of  the  non-linear  water  wave-wave  inter- 
actions. In  addition,  accurate  experimental  data  is  limited.  There- 
fore, a complete  theory  covering  all  the  stages  of  wind  wave  generation 
processes  has  not  yet  been  established. 

Ideally,  a comprehensive  theory  of  wind  wave  generation  should 
give  adequate  predictions  during  the  various  stages  of  wave  growth  for 
a given  fetch,  wind  speed,  and  direction.  In  other  words,  a complete 
theory  should  be  able  to  determine  the  state  of  the  sea  for  a given 
wind  field.  The  dynamical  wave  theory  is  concerned  with  the  local 
interactions  of  the  wave  field  in  a coupled  ocean-atmosphere  system. 
That  the  wind-wave  problem  is  similar  to  that  of  an  instability  between 
fluids  with  a density  discontinuity  and  moving  relative  to  each  other 
was  first  suggested  by  Helmholtz  (1868).  The  problem  was  treated  in 
more  detail  by  Thomson  (Lord  Kelvin,  1871).  The  well  known  "Kelvin- 
Helmholtz  Instability"  was  the  first  theory  to  explain  the  wind-wave 
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generation.  The  problem  was  further  studied  by  Jeffrey  (1925),  Wuest 
(1949),  Lock  (1954),  Eckart  (1953)  and  Ursell  (1956).  Recently,  the 
most  significant  contributions  to  solution  of  this  problem  were 
Phillips'  (1957)  and  Miles'  (1957)  theories  of  wave  generation  which 
considered  the  momentum  transfer  from  the  wind  field  into  the  waves. 
Phillips  (1960)  and  Hasselmann  (1960,  1962,  1963)  made  further  contri- 
butions to  the  nonlinear  wave-wave  interaction  theory  which  considered 
the  energy  transfer  among  different  wave  components  in  a wave  spectrum. 

Phillips'  (1957)  linear  theory  assumed  that  the  water  motions  are 
inviscid,  irrotational,  small  amplitude,  two-dimensional,  and  the  exci- 
tations of  the  water  surface  are  not  affected  by  the  wave-induced 
perturbation  in  the  air.  This  so-called  "resonance  model"  was  con- 
structed by  considering  turbulent  wind  flow  over  a water  surface  which 
is  initially  flat.  The  turbulent  pressure  fluctuations  in  the  wind 
field  generate  oscillations  of  the  water  travelling  in  all  directions 
and  having  all  the  wave  numbers  appearing  in  the  air  pressure  spectrum. 
When  a water  wave's  speed  matches  with  the  speed  of  translation  of 
the  corresponding  air  pressure  fluctuation,  the  wave  will  resonate  with 
the  pressure  field  and  grow  as  long  as  the  pressure  fluctuations  retain 
their  phase  relative  to  the  wave.  This  model  results  in  a linear  growth 
of  the  wave  spectrum  with  time  in  the  initial  growth  stage  and  gives 
order  of  magnitude  agreement  with  the  observed  wave  growth.  The  model 
becomes  invalid  when  the  wave  amplitude  is  sufficiently  large,  since 
the  neglected  wave-induced  perturbation  in  the  air  should  be  taken 
into  account  for  the  later  stage  of  wave  growth. 

Miles  (1957)  assumed  an  incompressible,  inviscid  air  shear  flow 
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over  single  frequency  progressive  water  waves.  He  considered  the  feed- 
back of  the  wave-induced  velocities  and  pressure,  but  the  turbulence  in 
the  air  was  neglected  except  its  role  in  establishing  the  turbulent 
logrithmic  velocity  profile.  The  water  motion  is  assumed  to  be  invis- 
cid,  irrotational , two-dimensional  and  of  small  amplitude.  The 
governing  equations  for  the  wave-induced  perturbations  in  the  air  stream 
are  exactly  the  same  as  those  in  the  laminar  instability  theory  (Orr- 
Sommerfeld  equation).  The  theoretical  result  is  that  the  momentum 
transferred  to  the  waves  by  the  action  of  the  wave-induced  Reynolds 
stresses  working  against  the  mean  velocity  gradient  is  evaluated  inside 
the  critical  layer  where  the  mean  wind  velocity  equals  the  wave  propaga- 
tion speed.  The  theory  is  applicable  for  those  wave  components  whose 
critical  layer  lies  outside  the  viscous  sublayer  and  in  the  region  where 
the  gradient  and  the  curvature  of  the  mean  velocity  profile  and  the 
magnitude  of  the  wave-induced  vertical  velocity  are  large. 

In  recent  years.  Miles'  theory  was  compared  with  experimental  data 
on  wave-induced  pressure  and  wave  growth  rate  obtained  both  in  labora- 
tories and  in  the  field.  Shemdin  and  Hsu  (1967)  measured  the  wave-induced 
pressure  distribution  over  mechanically  generated  waves  in  the  labora- 
tory and  inferred  that  the  wave  growth  rate  was  larger  by  a factor  of 
about  two  than  that  predicted  by  Miles'  theory.  Bole  and  Hsu  (1969) 
measured  the  wave  height  of  mechanically  generated  waves  along  the  wind- 
water  channel  and  found  that  the  growth  rate  was  underpredicted  by 
Miles'  theory  by  about  one-order  of  magnitude.  Similar  results  were 
obtained  in  the  field  by  Snyder  and  Cox  (1966),  Barnett  and  Wilkerson 
(1967),  and  Dobson  (1969).  Based  on  the  existing  experimental  results 

-3- 





we  can  conclude  that  Miles'  theory  is  inadequate  for  the  prediction  of 
wind-wave  growth.  Consequently,  investigators  have  been  devoting  their 
efforts  to  investigate  the  neglected  effects  of  turbulence  as  well  as 
nonlinear  mechanisms  on  wave  growth. 

Three  nonlinear  mechanisms  were  proposed.  First,  Miles  (1967) 
improved  his  own  model  by  considering  the  turbulence  which  was  neglected 
in  his  laminar  model.  The  expected  result  is  that  the  turbulent  struc- 
ture over  the  wavy  boundary  will  be  affected  by  the  interaction  with 
wave-induced  perturbation  field,  and  the  oscillating  wave  associated 
turbulent  stresses  may  be  essential  to  the  processes  of  energy  transfer 
to  the  water  waves.  Good  turbulence  and  wave-induced  perturbation  data 
are  needed  for  modeling  the  relationship  between  these  two  fields. 
Furthermore,  in  order  to  solve  the  modified  Orr-Sommerfeld  equation  con- 
taining the  unknown  terms  associated  with  the  oscillating  wave-induced 
turbulent  stresses,  ad-hoc  closure  conditions  would  have  to  be  assumed. 
Efforts  were  made  by  Hussain  and  Reynolds  (1970),  Davis  (1970,  1972)  and 
Chao,  et  al.  (1976).  As  we  know,  this  is  one  of  the  basic  unsolved  pro- 
blems associated  with  the  process  of  the  generation  of  turbulent 
stresses.  Therefore,  more  effort  is  needed  before  we  truly  understand 
the  role  of  the  oscillating-wave-associated  turbulent  stresses. 

Second,  Longuet-Hlggins  (1969a)  proposed  a nonlinear  "Maser  me- 
chanism" based  on  the  commonly  observed  phenomenon:  long  waves  sweeping 
over  short  waves  and  causing  short  waves  to  break  on  the  forward  face  of 
the  long  wave  crests  and  to  give  up  their  energy  to  the  long  waves.  The 
energy  transfer  rate  is  proportional  to  the  orbital  velocity  of  the  long 
waves  since  the  orbital  motion  is  forward  at  the  forward  face  of  the 
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long  wave  crest,  the  short  wave  therefore  supplied  a positive  amount  of 
energy  to  the  long  waves.  However,  Hasselmann  (1971)  proved  that  this 
positive  energy  transfer  is  almost  exactly  balanced  by  the  loss  of  po- 
tential energy  arising  from  the  mass  transfer.  Analyzing  the  field 
data,  Barnett  (1971)  failed  to  identify  the  energy  transfer  caused  by 
this  mechanism,  he  concluded  that  "maser  mechanism"  is  not  of  first 
order  importance  to  the  energy  balance  of  the  wave  spectrum. 

Third,  Hasselmann  (1962,  1963)  proposed  the  nonlinear  wave-wave 
interaction  theory  which  shows  that  the  energy  can  be  transferred 
among  different  wave  components  in  a gravity-wave  spectrum.  He  em- 
ployed perturbation  technique  to  find  the  nonstationary  growth  in  the 
initial  energy  spectrum.  Appreciable  modification  of  the  initial  sea 
state  was  found  from  the  solution  of  the  fifth-order  perturbation 
equation.  The  rate  of  change  of  the  energy  spectrum  can  be  expressed 
by  an  integral  expression  in  terms  of  wave  spectral  density  and  a 
transfer  function. 

Figure  1.1  (Barnett  1971)  shows  the  comparison  of  energy  transfer 
calculated  from  Hasselmann's  nonlinear  wave-wave  interaction  theory  by 
using  an  empirical  fully  developed  Pierson-Moskowitz  spectrum  and  ob- 
served field  spectrum.  For  the  broader  Pierson-Moskowitz  spectrum,  the 
major  transfer  acts  to  sharpen  the  spectral  peak,  while  the  nonlinear 
transfer  in  the  observed  spectrum  acts  to  maintain  the  low  frequency 
steep  forward  face  of  the  spectrum  and  serves  as  an  influential  growth 
mechanism.  The  figures  show  an  important  fact,  namely,  the  energy 
transfer  is  always  from  the  high  frequency  portion  to  the  low  frequency 
portion  of  the  spectrum.  Based  on  JONSWAP  (Joint  North  Sea  Wave 


Project)  data,  Barnett  (1971)  found  that  from  30%  to  90%  of  the  observed 
wave  growth  is  accounted  for  by  the  nonlinear  wave-wave  interaction 
theory.  Therefore,  Barnett  suggests  that  the  nonlinear  wave-wave  inter- 
action mechanism  will  dominate  all  the  energy  transfer  soon  after  the 
wave  spectrum  generated  by  Phillips'  resonance  mechanism.  Because  it 
takes  excessive  computer  time  to  integrate  the  triple  integral  equation 
for  the  nonlinear  energy  transfer,  Barnett  (1966)  proposed  a parame- 
terized equation  to  approximate  Hasselmann's  theoretical  results  for 
practical  applications.  This  parametric  equation  is  applicable  for  both 
fully  and  partially  developed  wave  spectra.  Mitsuyasu  (1968b)  employed 
Barnett's  parametric  equation  in  his  investigation  of  nonlinear  energy 
transfer  in  a wave  decaying  region  and  indicated  the  importance  of  this 
mechanism. 

Barnett  (1968)  developed  a wave  prediction  model  based  on  the  radia- 
tive transfer  equation.  This  model  takes  into  account  energy  transfer 
by  both  resonance  and  instability  mechanisms,  the  nonlinear  wave-wave 
interaction  mechanism,  and  the  wave-breaking  dissipative  mechanism. 
Although  it  was  only  a first  attempt  to  establish  a framework  for  ra- 
tional wave  prediction,  reasonable  results  were  obtained. 

It  will  be  shown  later  that  successful  prediction  of  a wave  spec- 
trum evolution  depends  not  only  on  the  source  function  (associa- 

ted with  the  nonlinear  wave-wave  interaction) , but  also  on  the  source 
function  S^n  (turbulent  air-input  to  waves)  and  S^s  (wave  dissipa- 
tion due  to  white-capping).  Hasselmann(1974)  proposed  a white-capping 
dissipative  model  which  could  be  useful  for  a wave  energy  balance 
scheme.  If  the  appropriate  air-input  model  could  be  constructed,  then 


we  should  be  able  to  formulate  a more  complete  wave  prediction  model 


based  on  the  radiative  transfer  equation  with  source  terms  of  air-input, 
nonlinear  wave-wave  interaction,  and  white-capping  dissipation. 

One  of  the  most  important  questions,  which  must  be  answered,  in 
the  wave  generation  process  is  "What  fraction  of  the  momentum  available 
in  the  air  can  be  transferred  directly  into  the  wave  field?"  Stewart 
(1961)  reanalyzed  some  older  wave  data  measured  at  known  fetches  and 
wind  speeds  and  came  to  the  conclusion  that  the  lower  limit  is  about 
0.2  and  the  upper  limit  is  unity.  Based  on  Miles'  theory,  Phillips 
(1966)  estimated  this  fraction  for  waves  with  c/uA  > 5 , where  c 
is  the  wave  speed  and  u^  is  the  friction  velocity,  and  found  the 
fraction  to  be  less  than  or  equal  to  0.1.  Such  a small  fraction  of 
momentum  flux  to  the  wave  field  must  be  associated  with  an  underestima- 
tion of  the  energy  transfer  in  Miles  theory.  Wu  Jin  (1968)  measured 
the  wind  generated  waves  in  the  laboratory  and  his  results  were  in  the 
range  from  0.1  to  0.4  for  cp/u*  < 7.5,  where  Cp  is  the  wave  speed  of 
the  peak  frequency  of  the  wave  spectrum.  Based  on  this  field  measure- 
ments, Dobson's  results  are  in  the  range  from  0.6  to  1.6  for 
9.6  < Cp/u*  < 49  . The  average  fraction  is  1.1,  which  is  greater  than 
1 and  seems  unreasonable.  One  of  his  most  reliable  sets  of  data  yield 
a value  of  0.8.  The  momentum  in  the  air  was  estimated  from  the  empi- 
rical  formula,  pfl  c^  u^  (where  pa  is  air  density,  the  drag  coeffi- 
cient cj  equals  0.0012,  and  u^  is  the  mean  velocity  above  five 
meters  of  the  water  surface).  The  uncertainty  in  drag  coefficient  may 
be  a main  source  contributing  to  the  large  apparent  error.  Because  of 
differences  in  the  range  of  Cp/uA  values,  comparisons  between  Dobson's 
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and  Wu's  experimental  data  can  not  be  made.  Barnett  (1971)  analyzed 
the  field  data  by  computing  the  atmosphere  external  source  function  in 
the  energy  balance  equation  and  found  that  the  average  values  of  the 
momentum  fraction  varied  from  0.1  to  0.2.  There  are  close  to  the 
Phillips'  prediction  and  Wu's  data,  but  deviate  considerably  from 
Dobson's  results.  The  recent  JONSWAP  results  presented  by  Hasselmann, 
et  al.  (1973)  showed  that  the  fractions  are  at  least  0.8  for  x - 10-  , 
and  0.2  for  10 ^ <.  x < 104  , where  x is  nondimensional  fetch 
(x  = xg/u‘1Q  , x is  fetch,  g is  the  gravitational  acceleration  and 

u is  the  mean  velocity  at  ten  meters  above  the  water  surface) . 

10 

Therefore,  these  fraction  could  be  fetch  and  wind  speed  dependent. 
Careful  and  direct  measurement  of  this  fraction  must  be  done  in  order 
to  supply  correct  information  to  establish  a complete  wave  generation 
theory. 

1.2  Objectives 

Based  upon  the  previous  discussion,  we  conclude  that  Miles'  model 
alone  can  not  predict  wave  growth  successfully.  The  relative  importance 
of  the  nonlinear  wave-wave  interaction  mechanism  should  be  investigated, 
especially  as  the  basis  for  a practical  wave  prediction  scheme.  There 
have  been  no  systematic  experiments  designed  to  investigate  this 
mechanism  in  the  laboratory.  Although  there  are  some  preliminary  field 
data  on  nonlinear  transfer,  the  limited  accuracy  of  the  field  data  makes 
the  interpretation  of  the  experimental  results  difficult,  leading  to 
some  reservation  about  the  results.  It  seems  to  be  expedient  to  verify 
presence  of  the  non-linear  wave-wave  interaction  process,  to  obtain  a 
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quantitative  determination  of  the  momentum  and  energy  transfer  from 
wind  to  waves,  and  to  examine  Hasselmann's  proposed  dissipation 
mechanism  due  to  white-capping  of  a gravity  wind  wave  spectrum  in  a 
laboratory  under  controlled  conditions.  Therefore,  a complete  experi- 
mental program  was  proposed  and  carried  out.  The  main  objectives  are 
as  follows: 

(1)  To  obtain  the  simultaneous  experimental  data  of  wave  height, 
fluctuating  pressure  and  velocity  components  in  close  proximity  to  a 
wave-perturbed  air-water  interface  in  order  to  evaluate  the  momentum 
and  energy  transfer  from  wind  to  waves. 

(2)  To  study  and  verify  Hasselmann's  nonlinear  wave-wave  inter- 
action mechanism  experimentally,  and  to  verify  Barnett's  approximate 
parameterized  equation  for  practical  applications  in  calculating  energy 
transfer  due  to  nonlinear  wave-wave  interactions. 

(3)  To  study  and  verify  Hasselmann's  white-capping  dissipative 


mechanism  experimentally. 


2 . THEORETICAL  BACKGROUND 


j 
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The  energy  balance  equation  of  a gravity  wind  wave  field  during 
its  generation  processes  is  described.  The  relevant  source  functions 
associated  with  the  wave  generation  are  also  presented. 

2.1  Energy  Balance  Equation  of  Wind  Wave  Spectrum 

The  energy  balance  equation  (radiative  transfer  equation)  for  deep 
water  gravity  wave  was  proposed  independently  by  Gelci,  et  al.  (1956), 
Hasselmann  (1960)  and  Groves  and  Melcer  (1961) . This  equation  describes 
the  propagation  and  the  processes  of  generation  and  dissipation  of  a 
wind  wave  spectrum  and  can  be  represented  by  [The  more  complete  form 
formulated  by  Hasselmann  (1968)  is  presented  in  Appendix  B.] 

^ = |^(f,e,x,t)  + v(f ,0)  vF(f,e,x,t)  = S'(f,e,x,t) 

” S^n(f»0,x,t)  + S^1(f,0,x,t)  + S^s(f,0,x,t) 

+ S^.e(f,0,x,t)  (2.1) 

where  F is  the  two-dimensional  (frequency  f and  angle  0 ) energy 
spectrum  defined  by  the  ensemble  average  (denoted  by  the  bracket  < > ) 
of  the  wave  height  n , i.e., 

j'"  F(f,0)  dfd0  - r E(f ) df  - <n2>  (2.2) 

—n  o "o 

and  V(f,0)  is  the  group  velocity  of  the  wave  with  frequency  f flnd  pro- 
pagating at  an  angle  0 measured  clockwiae  from  the  wind  (x-)  direction. 
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In  equation  (2.1),  s'  is  the  net  rate  of  energy  transfer  due  to  all 
the  interaction  processes  existing  in  the  fields  of  wind,  wave,  current, 
etc.  is  the  rate  of  energy  transfer  from  wind  to  waves,  which 

includes  the  linear  mechanisms  proposed  by  Phillips  (1957)  and  Miles 
(1957),  and  nonlinear  interactions  between  the  wave-induced  field  and 
background  turbulent  field.  is  the  rate  of  energy  transfer  re- 

sulting from  Hasselmann's  (1962,  1963)  nonlinear  wave-wave  interactions. 

is  the  rate  of  energy  dissipation  caused  by  white-capping;  the 
theoretical  model  of  this  transfer  was  proposed  by  Hasselmann  (1974). 

S'  is  the  rate  of  energy  transfer  between  wind-mean  current  and  wave- 
mean  current  fields,  etc. 

Because  of  the  extremely  limited  variation  of  energy  spectra  with 
respect  to  0 in  laboratory  channels,  one-dimensional  (frequency  f) 
energy  spectra  were  measured  and  used  in  the  equation.  Consequently, 
equation  (2.1)  should  be  simplified  to  one-dimensional  form  for  this 
study.  It  is  commonly  assumed  that  under  the  assumption  of  the  in- 
dependence of  frequency,  f and  direction,  9,  F(f,6)  can  be  approxi- 

mated by  the  product  of  the  one-dimensional  energy  spectrum  E(f)  and 
the  normalized  directional  spreading  factor  D'  (0),  expressed  as 

F(f ,0)  = E(f ) D'(0)  (2.3) 


where  E(f)  is  defined  by  equation  (2.2)  and  D'  (9)  is  chosen  (Hassel- 
mann, 1963b;  Mitsuyasu,  1968b)  for  the  narrow  directional  distribution 
of  a wind-wave  channel  as 
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D ' (0 ) = (8/3ir)  Cos40 
= 0 

and  j D (0)d0  ■ 1 


for  | 0 | <_  tt 

for  1 0 1 > it 

(2.4) 

Substituting  equations  (2.3)  and  (2.4)  into  (2.1)  and  integrating  with 
respect  to  0 from  -it/ 2 to  tt/2  , by  assuming  the  absence  of  backward 
scatter  in  the  wind-wave  channel,  we  can  reduce  equation  (2.1)  to  a one- 
dimensional  equation  (frequency  dependent  only).  After  the  integration, 
equation  (2.1)  can  be  further  simplified,  in  the  x-direction  along  the 
channel  and  for  a stationary  deep  water  gravity  wave  spectrum,  to 


Pw  * V£)  ^ " S(f)  " Sin(f)  + Snl(f)  + Sds(f>  (2‘5) 

— „1f/2  2 

where  V (f)  = J V(f)  Cos0  F(f,0)  d0  / E(f)  = 32g  / (45it  f)  is  the 

x — tt/2 

mean  group  velocity  in  x-direction,  V(f)  = g /4irf  is  the  group  velo- 
city of  the  wave  with  frequency  f and  is  equal  to  half  of  the  phase 

speed  in  deep  water  and  p is  the  water  density.  Because  of  the  small 

w 

drift-current  in  the  channel  and  weak  interaction  between  wave  and 
current,  Sfe  in  equation  (2.1)  is  neglected.  S^(f),  Sn^(f)  and  Sjs(f) 
in  equation  (2.5)  represent  the  first  order  processes  of  the  energy 
transfer  associated  with  a wave  spectrum.  The  detailed  description  of 
these  transfer  processes  is  given  in  the  following  sections. 


j-L 


•—  • »- 
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2.2  Energy  Transfer  from  Wind  to  Waves 


The  energy  transfer  from  wind  to  waves  in  a turbulent  developing 
boundary  layer  is  derived  in  Appendix  D.  The  work  done  by  the  normal 
pressure  in  quadrature  with  the  wave  height  (Stewart  1961)  is 

Sln">  ' -■><'>  TT1  <2'6) 

where  p(t)  is  the  pressure  at  the  interface,  3n(t)  / 3t  is  the 
vertical  velocity  of  the  water  surface.  [The  negative  sign  is  associa- 
ted with  a defined  positive  upward  cartesian  co-ordinate  system.] 
Because  of  the  limitation  of  the  present  laboratory  instrumentation  and 
the  difficulties  of  following  the  whole  wave  spectrum,  it  is  exceeding 
difficult  to  measure  the  pressure  at  the  interface.  The  alternative  is 
to  use  a fixed-frame  of  reference  and  to  measure  the  pressure  as  close 
to  water  surface  as  possible. 

The  energy  transfer  spectrum  can  be  obtained  by  taking  the  Fourier 
Transform  of  equation  (2.6),  viz., 

Sln(f)  - 2»f  Qu  (f)  (2.7) 

where  Qup^(f)  is  the  quadrature  spectrum  of  pressure  and  wave  height. 
Since  the  ratio  of  energy  to  momentum  for  a deep  water  gravity  wave  is 
equal  to  its  phase  speed  c,  Stewart  (1961)  suggested  the  spectrum  of 
momentum  transfer  can  be  expressed  as  (see  also  Dobson,  1971) 

Tw(f)  “ Sin(f)  ' c “ 2lrf  Sln(f)  ' g (2.8) 
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where  c = g / 2^f  is  the  phase  speed  of  the  deep  water  gravity  wave. 
Based  on  equations  (2.7)  and  (2.8),  we  can  evaluate  the  energy  and 
momentum  transfer  from  wind  to  waves  due  to  normal  pressure  acting  on 
water  surface. 

2.3  Energy  Transfer  of  Nonlinear  Wave-Wave  Interactions 

2.3.1  Theoretical  Expression  of  Energy  Transfer  of  Nonlinear  Wave- 

Wave  Interactions 

Hasselmann  (1962)  proposed  the  conservative  nonlinear  wave- 
wave  interactions  mechanism  for  the  gravity  wind  wave  spectrum.  As  a 
result  of  this  interaction,  energy  can  be  transferred  and  redistributed 
among  different  wave  components.  The  theoretical  derivation  was  based 
on  the  assumptions  of  inviscid  fluid,  two-dimensional  and  irrotational 
water  motion.  The  velocity  potential  $(x,t)  and  free  surface  n(x,t) 
were  used  as  the  perturbation  parameters.  Hence  they  were  expanded 
in  a series  and  substituted  in  the  governing  Laplace  equation,  kinema- 
tical  and  dynamical  boundary  conditions,  as  well  as  in  initial  condi- 
tions to  generate  higher  order  perturbation  equations.  The  aim  of  his 
investigation  was  to  determine  the  change  (nonstationary  growth)  in  the 
initial  energy  spectrum  brought  about  by  the  nonsteady  higher-order 
perturbations.  On  using  a fifth  order  perturbation,  the  nonstationary 
solution  was  found  to  grow  continually  with  time,  resulting  in  an 
appreciable  modification  of  the  initial  sea  state.  The  rate  of  change 
of  the  energy  spectrum  was  determined  by  the  following  equations 
(Hasselmann  1963). 


oo  oo  7T 


fr[f(k4*e4)]  mS0S0  I.v  ^ w4  f<*2’0  2 ' f ^3.6^3) 


+ u>3  f(k1>e1)  f(k2,0vJ;)  f(k4,e4) 


- MjCkj/l^)  f(k3,e^;)  f(k4,e4)  f^.e^ 


“ifk]./^)  f Oc3,e(^))  f(k4,e4)  f(k2,e(^)]}  de1dk1dk2 


where  k and  w are  the  wave  number  and  frequency,  respectively. 
The  resonance  conditions  are 


ki  + £2  = £3  + fc4 


W.  + 0J»  *=  0)-  + 0). 
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The  dispersion  relation  for  a deep  water  wave  is 


iDj  ■ 2rrf j - (gkj ) 


The  equations  of  resonance  conditions  and  dispersion  relation  can  be 
reduced  to  the  single  equation  (Hasselmann,  1963b) 


i 


(2.11) 


2 CosB  |k2|  |kx  - k4|  - (/k^  + - ^)4  - k22  - (kx  - k4)2 

-V  -4-  -f 

with  6 being  the  angle  between  k2  and  (k^  - k^).  Since  the 

-4  -4 

interaction  surface  is  three  dimensional  in  (k^,k2)-space,  the  trans- 
formation to  the  surface  integral  can  be  affected  simply  by  introducing 

-4 

polar  coordinates,  i.e.  k^  - (k^,©^).  In  examining  equation  (2.9) 
for  fixed  k^  and  0^,  and  the  given  values  of  the  integral  variables, 
k^,  k2>  and  0^,  it  is  seen  that  ©2,  k^  and  0^  are  the  remaining 
undetermined  variables  for  completing  the  integration.  8 can  be 

determined  from  equation  (2.11)  by  giving  the  values  of  k^,  0^,  k^ , 

-*■  -*■  -*■ 

©4  and  k2>  Since  B is  the  angle  between  k2  and  (k^  - k^) , 
can  be  determined  by  providing  the  values  of  k2 , k^,  k^,  0^,  ©4 
and  B.  Finally,  k^  and  0^  can  be  calculated  by  equation  (2.10) 
for  all  determined  values  of  k^,  0^,  k2>  ©2,  k^,  and  0^.  ©2^\ 

0^^  and  are  correspondent  to  the  values  of  +B  and 

-B  respectively.  The  spectral  density  in  polar  co-ordinates  is 


f (k,0)  - P if  F(f ,0)  / 2k 
w 


and  the  function  T can  be  expressed  as 


T . 3gP  2 __3 

2pw  V2W3W4  k2  |kx  - kj  I Sine  | g2 


for  CosB  < 1 


for  CosB  > 1 
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where  D - D+  ■*  £ is  the  transfer  coefficient.  The  full  expres- 

l * *^2  * 

sions  for  D were  given  by  equations  (4.9)  and  (4.10)  in  Hasselmann 
(1962)  . Based  on  the  described  relationship  of  all  6 and  k and 
the  defined  functions  of  f(k,0)  and  T,  we  can  calculate  the  rate  of 
nonlinear  energy  transfer  due  to  water  wave-wave  interaction  by 
equation  (2.9). 

The  results  of  nonlinear  energy  transfer  (equation(2.9))  were 
interpreted  by  Hasselmann  (1962)  in  terms  of  quadruple  interactions 
between  three  "active"  wave  components  which  determine  the  interaction 
rate,  and  a "passive"  fourth  component  which  receives  energy  from  the 
first  three  components  but  has  no  direct  influence  on  the  interaction. 
The  numerical  calculations  of  a fully  developed  Neumann  spectrum  were 
performed  by  Hasselmann  (1963b).  The  results  showed  that  the  transfer 
process  tends  to  reduce  sharp  peaks  in  the  spectrum  and  redistribute 
the  energy  more  uniformly  over  all  wave  components  by  transferring  the 
energy  from  the  middle  frequency  range  into  the  lower  and  higher  fre- 
quency ranges. 

2.3.2  Parametric  Approximation  of  Energy  Transfer  of  Nonlinear 
Wave-Wave  Interactions 

The  evaluation  of  the  nonlinear  wave-wave  interactions 
involves  complicated  numerical  computations  with  equation  (2.9),  re- 
quiring six  hours  computer  time  (Hasselmann  1963b)  on  a CDC  1604. 

Even  the  newly  developed  program  by  JONSWAP  still  takes  20  minutes  of 
computer  time  on  a CDC  6600.  Thus,  it  seems  to  be  uneconomical  and 
impractical  for  field  application.  In  addition,  the  new  program 
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provided  by  Dr.  Sell  can  not  be  directly  used  for  this  study  because  of 
the  differences  in  frequency  scales.  Efforts  were  made  to  change  to 
frequency  scales  consistent  with  this  study,  but  unfortunately,  this 
did  not  succeed.  Consequently,  Barnett's  (1966)  parametric  equations 
to  approximate  Hasselmann's  theoretical  results  were  used  in  this  study. 

The  parameters  to  characterize  the  wave  field  are  defined  as 

TT  °° 

Energy:  e =J  J F(f,0)dfd0 

-it  o 

l/2 

Mean  frequency:  f = [—  f J F(f,0)  f2  dfd0] 

o e * -tj  " 0 
I TT  00 

Mean  direction:  0 = — f f F(f,0)  0 dfd0 

° e -it  0 


The  parametric  equation  of  energy  transfer  associated  with  non- 
linear wave-wave  interactions  is 

s'  (f  ,0)  = [T'(f,0)  - t’  (f,0)  F(f ,0) ] pc  (2.14) 

ni  w 

where  T and  t'  are  integral  functions  of  F.  This  approximation 
was  considered  to  be  the  sum  of  two  processes  as  originally  demon- 
strated by  Hasselmann  (1962).  r'  was  related  to  the  passive  components 
which  receive  energy  from  other  components,  while  t'  F was  related  to 
the  active  components  which  transfer  energy  to  other  components.  Using 
dimensional  analysis  and  the  qualitative  features  of  Hasselmann's 
theory,  Barnett  proposed  the  following  forms  of  r'  and  t'  for  the 
case  of  fully  developed  sea: 
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(f.9)  ■ Co84( 


9.  9 ) 

o'  v f ' 


• exp [-4(1  - -p)2  + 0 . 1 Op.) 5 ] 


for  f > 0.42f  and  |e  - 0 | < ir/2 


■ 0 otherwise 


(1  + 16  | Cos(0  - eo)|)  fo7  (f  - 0.53f„)3 


for  f > 0.53f„ 


” 0 otherwise 


The  one-dimensional  form  of  the  nonlinear  energy  transfer  equation 
can  be  derived  by  substituting  equations  (2.4),  (2.15)  and  (2.16)  Into 
equation  (2.14)  and  Integrating  it  with  respect  to  6 from  — tt/2  to 
n/2.  The  resulting  form  can  be  expressed  as 


Snl(f)  = tF(f)  " T(f)  E(f)1  pw8 


rm  . ^•^•108e3f»8  3*  ,f  - 0.42fo,3 

K ' „4  8 1 f ; 


• exp [-4(1  - |^)2  + 0.1(|V] 


for  f > 0.42fo  and  |e  - 60|  < tt/2 


0 otherwise 


(2.18) 


tm 


and 


T(f)  . ■7.-.5\lp7.e2  (1  +I6^8)  ^7  (f  _ 0>53fo)3 


for  f > 0.53fg 


■ 0 otherwise 


(2.19) 


The  results  of  nonlinear  energy  transfer  calculated  with  the  para- 
meterized equation  (2.17)  compared  well  with  those  of  Hasselmann ' s 
theoretical  equation  for  the  fully  developed  and  partially  developed 
seas.  The  comparison  for  a typical  fully  developed  sea  is  shown  in 
Figure  2.1. 

2.4  Energy  Dissipation  of  White  Capping 

Because  dissipative  processes  such  as  molecular  viscosity  and 
turbulence  (Phillips  1959,  Hasselmann  1968)  are  inadequate  to  account 
for  the  observed  energy  removal  from  wind  waves,  the  white  capping  is 
generally  believed  to  be  the  dominant  dissipative  mechanism  in  a wave 
field  at  moderate  and  higher  wind  speeds.  The  effect  of  white  capping 
on  the  spectra]  energy  distribution  was  investigated  by  Hasselmann 
(1974).  The  problem  was  treated  by  expressing  the  white-cap  inter- 
actions in  terms  of  an  equivalent  ensemble  of  random  pressure  pulses. 

It  is  first  shown  that  the  source  functions  S.  for  any  non-expansible 

ds 

Interaction  process  which  is  strongly  nonlinear  locally  but  still 
weak-in-the-mean  is  quasi-linear  with  respect  to  the  wave  spectrum, i .e . , 

s”  (k)  - -y'b'oo  (2.20) 

as 
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where  Y is  the  damping  factor  which  depends  on  the  entire  wave 
spectrum;  y'  is  found  to  be  y'  “ d (2irf)2,  under  the  assumptions  that 
the  space-time  scales  of  a white  cap  are  small  compared  with  the  wave- 
heights  and  periods  of  the  waves.  Here,  d is  a constant  for  a given 
wave  state  and  can  be  determined  indirectly  from  the  energy  and  momentum 
balance  in  the  high-frequency  equilibrium  range.  Finally,  the  dissipa- 
tion spectra  of  ocean  waves  due  to  white  capping  can  be  expressed  in 
the  frequency  domain  by  assuming  that  d is  independent  of  the 
direction  6 ; the  spectra  is  then 

S.  (f)  - -d  (2irf ) 2 E(f ) p g (2.21) 

as  w 

We  can  measure  the  net  rate  of  energy  transfer  S (equation  (2.5)) 
and  the  rate  of  energy  transfer  from  wind  to  waves  (equation (2. 7) ) . 

Based  on  equation  (2.5),  we  can  subtract  from  S to  obtain 

experimentally  the  sum  S ^ + S^.  It  is  clear  that  in  order  to  sepa- 
rate the  sum  of  the  experimentally  determined  S ^ , either  S^s 

or  S , would  have  to  be  assumed.  Based  on  the  measured  wind  wave 
nl 

spectrum  in  the  channel,  equations  (2.17)  and  (2.21)  were  used  alter- 
nately in  calculating  the  Sni  and  Sjg  , respectively. 

2 . 5 Mean  Momentum  Transfer  from  Wind  to  Waves 

The  momentum  transfer  across  the  wind-wave  surface  in  a turbulent 
developing  boundary  layer  is  derived  in  Appendix  C.  By  neglecting  the 
relatively  small  viscous  terms,  the  final  results  of  momentum  transfer- 
red into  current  and  waves,  in  the  direction  of  wave  propagation,  can 


be  expressed  as 


- P; 


urvr 


(2.22) 


Pa 


u'2  |a 
3x 


(2.23) 


where  Mc  and  are  momenta  transferred  to  the  current  and  waves 

respectively.  From  these,  we  find  the  momentum  transferred  into  waves 
is  not  due  only  to  the  normal  air  pressure  acting  on  water  surface  but 
also  due  to  the  oscillating  turbulent  Reynolds  stresses  (wave  associated 
turbulent  stresses)  in  the  air.  Finally,  the  ratio  of  the  momentum 
transferred  into  the  waves  (or  supported  by  the  wave  form)  over  the 
total  momentum  transferred  across  the  interface  in  the  direction  of 
the  wave  propagation  can  be  estimated  by  the  following  equation 


Y - M„  / (Mc  + Mw) 


(2.24) 


* j 

*♦  V 
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3 . EXPERIMENTAL  APPARATUS 


Based  on  the  objectives  and  theoretical  background,  the  required 
data  are  water  wave  height,  fluctuating  pressure  and  velocities  in  the 
wind  stream  and  the  mean  ve' ■ jty  profiles  along  the  wind-wave  channel 
as  a function  of  fetch  for  a range  of  wind  speeds.  The  facilities  and 
instrumentation  relevant  to  this  experimental  study  are  presented  in 
the  following  sections. 

3 . 1 Channel  Description 

All  the  experiments  were  conducted  in  the  Stanford  Wind-Wave 
Channel  (Hsu,  1965)  in  the  Hydraulics  Laboratory  at  Stanford  University. 
A schematic  is  shown  in  Figure  3.1.  The  channel  is  37.7  m long  with  a 
rectangular  cross  section  which  is  2 m high  and  0.98  m wide.  The  glass- 
walled  test  section  is  approximately  23  m long.  Baskets  of  stainless 
steel  machine  turnings  are  located  at  the  end  of  the  channel  and  serve 
as  a beach  to  absorb  wave  energy.  An  airfoil-bladed  fan  located  at  the 
channel  exit  draws  air  through  the  test  section.  A honeycomb  is  placed 
at  the  end  of  the  channel  to  suppress  secondary  flows  caused  by  the 
centrifugal  action  of  the  fan.  The  air  inlet  contains  a curved  section 
with  a series  of  turning  vanes,  a 5.1  cm  thick  honeycomb  with  a 0.69  cm 
hexagonal  mesh,  and  a group  of  small  size  screens  to  produce  a relative- 
ly uniform  flow.  A 1.5  m long  flat  plate  was  mounted  at  the  base  of  the 
last  screen  to  provide  a smooth  transition  for  incoming  air  to  the  test 
section.  The  depth  of  water  was  approximately  one  meter  and  was 
adjusted  until  it  coincided  with  elevation  of  the  plate;  thus,  excessive 
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inlet  disturbances  were  eliminated  and  the  fetch  could  be  determined 


precisely. 

3.2  Instrumentation 

3.2.1  Mean  Velocity 

The  mean  velocity  profiles  were  measured  with  a 0.238  cm 
o.d.  United  Sensor  pitot-static  tube  connected  to  a Pace  differential 
transducer  (model  P90D)  with  a maximum  pressure  differential  of  0.037 
psid.  The  transducer  served  as  one  leg  of  a Sanborn  650  Recorder 
inductive  bridge  which  converts  pressure  differential  to  electrical 
signals.  A Combustion  Institute  micromanometer  with  resolution  of 
+ 0.00127  cm  of  silicone  oil  (r.d.  0.82)  was  used  to  calibrate  the 
transducer.  A typical  calibration  curve  is  shown  in  Figure  3.2.  The 
variation  of  the  slope  of  calibration  curve  was  less  than  3%  over  the 
past  three  years.  The  transducer  is  extremely  stable  and  reliable. 

The  mean  velocity  was  calculated  using  the  expression 

U = / 2Ap  / Pa  (3.1) 

where  Ap  is  differential  pressure  and  pa  is  the  local  air  density. 

3.2.2  Wave  Height 

The  elevation  of  the  instantaneous  water  surface  was 
measured  by  a capacitance  wave-height  gauge.  The  wire  gauge  is  an 
Isonel-insulated  0.1  mm  diameter  copper  conductor.  Variation  of  water 
levels  produces  a corresponding  change  in  the  capacitance  of  the 


insulted  wire.  A complete  description  of  the  wave  gauge  can  be  found  in 
Colonell  (1°66) . 

Calibration  of  the  wave  gauge  was  carried  out  by  varying  the  depth 
of  the  submergence  of  the  wire  in  the  still  water.  The  wave  gauge  was 
mounted  on  an  accurate  vertical  positioning  device  located  on  an  instru- 
mentation carriage  in  the  wave  channel.  The  distance  x(fetch)  variation 
can  be  achieved  by  moving  the  carriage  along  the  channel.  The  accuracy 
of  measuring  the  vertical  submergence  of  the  wave  gauge  is  0.254  mm. 
Figure  3.3  shows  a typical  calibration  curve  of  wave  height. 

3.2.3  Velocity  Fluctuations 

Thermo-Systems,  Model  1010  constant  temperature  anemometer 
and  a platinum  cross-hot-film  probe  (X-array)  having  a diameter  of 
0.0508  mm  were  used  to  obtain  the  fluctuating  velocity  components 
measurements.  The  frequency  response  of  the  cross-hot-film  (X-array) 
i 40,000  Hz.  The  longitudinal  and  vertical  turbulent  velocity  compo- 
nents were  determined  by  the  sum  or  difference  of  the  X-array  signals. 

The  hot  films  were  calibrated  against  the  pitot-static  tube  (Sec- 
tion 3.2.1)  in  the  free  stream  of  the  wind-wave  channel.  The  voltage 
output  E from  the  anemometers  and  the  effective  velocity  U ^ normal 
to  the  hot  film  can  be  represented  by 

E2  * A + B Un  (3.2) 

eff 

where  A,  B and  n are  constants  to  be  determined  from  the  calibra- 
The  effective  velocity  can  be  determined  by  resolving  the 


tion. 


measured  free  stream  velocity  normal  to  each  hot-film.  The  least 


tv.*. 


square  method  was  employed  to  determine  these  constants  for  each  film 

by  providing  the  data  of  E and  Ue^  for  ten  wind  speeds.  Although 

the  calibrations  were  found  to  be  quite  repeatable,  a check  of  the  wire 

calibration  was  carried  out  before  and  after  each  run.  A typical  cali- 
2 

bration  curve  of  E versus  U ,,  is  shown  in  Figure  3.4. 

ett 

3.2.4  Pressure  Fluctuation 

The  pressure  signal  was  measured  by  a crystal  pressure 
transducer  manufactured  by  the  Cornell  Aeronautical  Laboratory.  The 
assembly  consists  of  a thin  disc-shaped  pressure  probe,  connecting  tube, 
and  phenolic  disc-shaped  housing  which  contains  the  crystal  transducer 
sensor  and  high  impedance  converter.  The  output  from  the  high  impedance 
converter  was  connected  to  a 1000  gain  box  to  complete  the  pressure 
measuring  system.  The  probe  is  a 0.159  cm  thick,  0.635  cm  in  diameter 
streamlined  thin  disc  with  a 0.033  cm  diameter  piezometric  hole  at  the 
center  of  disc  on  each  side.  A small  stainless-steel  tube  of  the  same 
diameter  as  the  piezometric  hole  was  used  to  connect  the  probe  and  the 
transducer  sensing  element.  The  small  size  of  the  connecting  tubing  is 
necessary  to  minimize  the  inertia  effect  of  the  air  mass  in  the  system. 

A larger  (0.159  cm  diameter)  outer  tube  was  used  to  protect  the  inner 
tube  and  also  provide  better  supporting  strength.  The  detailed  struc- 
ture of  this  transducer  is  shown  in  Figure  3.5. 

The  crystal  pressure  transducer  has  a linear  input-output  response 
and  a sensitivity  of  2.0  volts  per  psi  with  a noise  level  of  50  micro- 
volts peak-to-peak . The  output  of  this  transducer  is  very  sensitive  to 
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the  effect  of  dirt  and  humidity  so  special  care  was  taken  to  maintain 
its  normal  operation. 

The  pressure  transducer  was  calibrated  dynamically  against  the 
Pace  transducer,  while  the  Pace  transducer  was  calibrated  statically  by 
a combust  micromanometer  (Section  3.2.1).  The  dynamic  calibration  was 
performed  in  a pressure  chamber  1.32  m long  with  a 12.7  cm  inner  diame- 
ter. A sinusoidal  pressure  signal  was  generated  by  the  reciprocal 
motion  of  a small  piston-cylinder  assembly  mounted  on  the  end  of  the 
pressure  chamber  and  connected  to  a rotating  wheel  driven  by  a D.  C. 
servo  motor  through  a D.  C.  function  generator.  During  calibration 
both  the  pressure  sensing  probes  connected  to  the  crystal  sensor  and  to 
the  reference  Pace  transducer  were  placed  at  the  end  of  the  pressure 
chamber  opposite  to  the  piston-cylinder  assembly.  Three  different 
amplitudes  of  the  sinusoidal  input  signal  can  be  obtained  by  changing 
the  stroke  of  the  piston.  A typical  calibration  curve  of  the  crystal 
transducer  output  versus  the  Pace  transducer  output  is  shown  in  Figure 
3.6.  The  frequency  response  of  the  crystal  pressure  transducer  is 
constant  over  the  range  of  1 to  20  Hz  which  was  satisfactory  for  this 
present  study.  Figure  3.7  shows  the  frequency  responses  at  the  three 
different  amplitudes.  The  magnitude  of  the  transfer  function 
( Ih C f ) I)  associated  with  the  pressure  sensing  system  was  defined  in 
Section  (4.2.3).  The  phase  shift  of  the  pressure  sensing  system 
(including  the  1000  gain  box)  was  calibrated  against  a pulse  signal 
triggered  by  the  light  passing  through  a hole  in  the  rotating  wheel  and 
to  a photo-electric  cell  located  at  the  other  side  of  the  wheel.  Both 
the  pressure  and  pulse  signals  were  recorded  simultaneously  on  a Sanborn 


recording  oscillograph  to  deduce  the  phase  shift  of  the  pressure  sensing 
system.  Figure  3.8  shows  the  phase  shift  induced  by  the  crystal  trans- 
ducer, 1000  gain  box  and  high  impedance  converter  as  a lag  of  the 
sensing  system  behind  the  actual  pressure  wave. 

3.2.5  Sensor  Arrangement 

The  wave-height  gauge,  the  hot  film  and  the  crystal  pressure 
sensor  were  installed  in  the  channel  at  the  same  fetch.  The  wave- 
height  gauge  was  installed  at  the  center  line  of  the  channel.  The  hot 
film  and  the  pressure  sensor  were  located  1.5  cm  to  each  side  of  the 
wave-height  gauge  and  mounted  on  a vertical  positioning  device  which 
can  be  adjusted  by  a control  panel  outside  the  channel. 

3.3  Traversing  Mechanism 

In  order  to  measure  the  mean  velocity  profiles  for  each  stations, 
the  pitot-static  tube  was  mounted  on  the  instrument  carriage.  The 
carriage  and  its  traversing  mechanism  could  be  moved  in  either  hori- 
zontal or  vertical  directions.  The  vertical  positioning  of  the  tra- 
verse was  remote  controlled  with  a resolution  of  0.25  mm.  The  details 
of  the  carriage  and  traversing  mechanism  are  given  by  Chambers,  et  al 
(1970). 

3.4  Data  Acquisition-Reduction  System 

A Hewlett-Packard  2100A  computer  system  was  used  for  data  acquisi- 
tion and  analysis.  This  real-time  data  acquisition-reduction  system 
includes  a computer,  analog-to-digital  converter,  magnetic  tape  drive. 


output  terminal,  disc  operating  system,  etc.  The  dynamic  range  of  the 
analog-to-digital  converter  is  + 10  volts.  Sixteen  active  input  chan- 
nels to  the  converter  are  available.  The  maximum  sampling  rate  is 
45,000  samples  per  second.  The  accuracy  of  the  system  is  + 0.01  volts 
with  a resolution  of  + 5 mv.  The  statistical  evaluation  of  this  system 
performance  was  made  by  Takeuchi  and  Mogel  (1974),  and  the  excellent 
agreement  with  theoretical  predictions  was  obtained. 

All  the  data  taken  by  this  system  were  digitized  and  stored  on 
magnetic  tape  for  later  analysis. 


4.  DATA  ACQUISITION  AND  REDUCTION  METHODS 


4 . 1 Data  Acquisition  Procedure 

In  this  study,  experimental  data  were  taken  at  five  stations  along 
the  channel  with  fetches  of  3.46  m,  6.51  m,  9.48  m,  12.61  m and  15.66  m 
and  three  wind  speeds  of  7.09,  8.01  and  8.88  m/sec.  The  channel  was 
filled  with  water  at  least  two  days  before  each  experiment  to  allow 
ample  time  for  deaeration.  Otherwise  the  dissolved  air  in  the  water 
adheres  to  the  wire  of  the  wave  height  gauge  and  affects  its  normal 
operation . 

The  normal  operating  procedure  was  to  set  the  desired  wind  speed 
and  run  the  channel  for  approximate  thirty  minutes  before  data  acquisi- 
tion to  insure  development  of  a stationary  wave  height  spectrum.  The 
highest  point  of  the  wave  was  then  determined  and  the  hot  film  and 
crystal  pressure  transducer  probe  was  set  0.508  cm  above  this  highest 
point.  The  fluctuating  pressures,  velocities,  and  wave  height  data 
were  simultaneously  taken  and  recorded  by  the  HP  2100A  data  acquisi- 
tion-reduction system.  All  the  signals  were  appropriately  amplified 
and  low  pass  filtered  with  the  3 db  point  at  20  Hz;  they  were  then 
sampled  every  0.025  seconds  for  10  minutes  and  stored  on  magnetic 
tape.  The  measurements  were  made  at  five  different  elevations 
(0.058  cm  apart  each)  above  the  highest  point  of  the  wave.  The  com- 
putation of  the  experimental  results  was  based  on  the  data  taken  at 
the  lowest  elevation,  i.e.,  0.508  cm  above  the  highest  wave.  The 
remaining  data  was  used  as  a check  to  reveal  any  abnormality  during 
the  data  taking  process.  However,  in  order  to  reveal  the  full 


frequency  range  of  the  turbulence  characteristics  in  the  air  stream, 
one  extra  set  of  velocity  fluctuation  data  was  taken  at  the  lowest 
elevation  with  a 5000  Hz  sampling  rate  and  2500  Hz  cutoff  frequency 
for  four  minutes.  There  were  six  runs  (five  elevations  and  one  extra 
velocity  at  lowest  elevation)  for  each  wind  speed.  After  finishing 
experiments  at  the  three  selected  wind  speeds,  the  instrumentation  was 
moved  to  another  station  until  data  taking  at  the  five  stations  was 
completed. 

The  mean  velocity  profiles  were  measured  independently  by  pitot- 
static  tube  after  completing  the  data  taking  discussed  previously. 

The  mean  velocity  was  obtained  by  averaging  the  readings  which  were 
taken  every  0.01  second  for  two  minutes.  The  experimental  conditions 
were,  of  course,  kept  the  same  as  those  used  previously. 

4.2  Data  Reduction  Procedure 

4.2.1  Auto  and  Cross  Spectra,  Coherence  and  Phase 

The  auto  and  cross  spectra  of  wave  height,  wave-induced 
velocities  and  pressure  were  computed  via  FFT  method  with  a cosine 
window  and  a 512  point  block  size  for  a 50  spectra  average.  The  auto 
spectrum  S^(f)  of  the  random  variable  x(t)  can  be  estimated  as 

Sx(f)  = <X*(f)X(f)>  (4.1) 

* 

where  X(f)  is  the  Fourier  Transform  of  x(t)  and  X (f)  is  the 
complex  conjugate  of  X(f).  The  bracket  represents  the  ensemble 
average  or  spectra  average,  which  is  the  averaging  process  over  the 


total  number  of  measurements  of  X (f)X(f)  at  each  frequency. 
The  cross-spectrum  of  the  variables  x(t)  and  y(t)  is 


S (f)  = <X  (f )Y(f ) > - Co  (f)  + 1 Qu  (f) 
xy  xy 


where  Co^^(f)  Is  the  cospectrum  and  Quxy(f)  is  the  quadrature 
spectrum. 

The  coherence  and  phase  spectra  can  be  defined  as 


2 2 1/2 
Rxy(f)  = ([Co  xy(f)  + Qu  xy  (f)  ] / [Sx(f)-Sy(f)]}i/‘i  (4.3) 


0 (f)  = tan  [Qu  (f)  / Co  (f)] 

xy  xy  xy 


(4.4) 


4.2.2  Auto  Spectrum  of  Wave-Induced  Velocity  Fluctuations 

The  measured  velocity  signal  x(t)  can  be  considered 
(Chang,  Plate  and  Hidy,  1971)  as  the  sum  of  the  mean  quantity  x,  the 
wave-induced  perturbation  x(t),  and  the  background  turbulence  x'(t). 
The  Fourier  Transform  of  x(t),  after  removal  of  the  mean  part  x,  can 
be  expressed  as 


X(f)  = X(f)  + X (f) 


(4.3) 


where  X(f)  is  the  quantity  correlated  with  the  wave,  and  X (f)  is 
the  random  component. 

The  auto  spectra  of  wave-induced  velocities  can  be  derived  under 
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the  assumption  of  the  statistical  independence  of  background  turbulence 
and  wave  height  n (Chang,  Plate  and  Hidy,  1971),  i.e.  <x'  (f)1f(f)>  =0, 
where  the  bracket  denotes  an  ensemble  average  and  TT(f)  is  the  Fourier 
Transform  of  n(t).  The  derivation  is  as  follows: 


Svn'Snv  = <X*JT>-<ir*X>  = (CO  + i Qu  ) (Co  - i Qu  ) 
XTi  nx  xn  xn  xn  xn 


2 2 

Co  + Qu 
xn  xn 


(4.6) 


f * 

By  appling  the  assumption  <X  1T>  = 0,  equation  (4.6)  can  be  written  as 
follows  after  we  take  a sufficient  number  of  spectrum  averages  (a  50 
spectra  average  was  used  for  this  data  reduction) : 


Sxn’Snx  = <X*TT>,<TT*X>  - <X*TT>-<TT*X> 


- <x  x>-<n  n>  - s-*s 

X T) 


(4.7) 


Equating  equations  (4.6)  and  (4.7),  one  finds  the  auto  spectra  of  the 
wave-induced  velocities  as 


S,(f)  » <X*(f)2(f)>  - (Co2  (f)  + Qu2  (f))  / S (f) 
x xn  xn  n 


(4.8') 


Another  technique  could  be  employed  for  this  particular  analysis 
and  the  same  result  obtained.  If  one  assumes  that  the  X(f)  response 
to  TT(f)  is  a linear  process,  i.e.,  X(f)  * X (f)  1T(f),  the  randomness 
in  X(f)  is  similar  to  that  in  11(f).  The  factor  [or  the  transfer 
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function]  x(f)  may  contain  information  of  phase  shift  and  amplitude 


attenuation  of  X(f)  with  respect  to  H(f).  If  one  multiplys  h'*(f) 
with  X(f)  in  the  assumed  linear  relationship,  one  obtains 

X(f)  TT*(f)  = X(f ) TT(f ) U*(f) 


Taking  the  ensemble  average  of  the  above  equation  produces 


<X(f)  T*(f)>  = X(f)  <TT(f)lT*(f)> 


Solving  for  X(f),  one  obtains 


. <X(f)fl  (f)> 
<ir(f)Ti*(f)> 


Consequently,  the  wave  related  signal  X(f)  can  be  expressed  as 


x(f)  = <x(f)l^(f)>  n(f) 
<11  ( f )tt  (f)> 


■ <X.Cff)n  .<f>>  ,1(f) 

<n(f)n  (f)> 


(4.9) 


t * 

since  <X  (f)TT  (f)>  = 0.  Based  upon  equation  (4.9),  the  averaged  power 


spectrum  of  a wave-induced  quantity  can  be  expressed  as 


-34- 


— - — 


■ 


s-(f)  - <X*(f)X(f)> 


/ * .T* 

/<x 

\ <^T*TI> 


<XT]*>-tl\ 

<rr*rr>  / 


* * / 
= <X  IT >•  <Xfl  > / 


<ir*rr> 


= (Co2  (f)  + Qu2  (f))  / S (f) 

Xfl  Xn  h 

This  expression  is  the  same  as  equation  (4.8),  which  reveals  that  the 
linear  assumption  is  reasonable  and  acceptable.  This  assumption  will 
be  employed  latter  for  analyzing  the  energy  and  momentum  transfer  from 
wind  to  waves  due  to  the  oscillating  wave-associated  turbulent 
stresses . 

4.2.3  Auto  Spectrum  of  Wave-Induced  Pressure  Fluctuation 

The  desired  pressure  signals  p(t)  are  modified  by  a 
transfer  function  h(t)  of  the  system  which  includes  the  crystal 
transducer,  high  impedance  converter  and  the  1000  gain  amplifier. 
Therefore,  the  measured  pressure  signals  Pm(t)  equal  the  convolu- 
tion of  p(t)  and  h(t),  viz., 

P (t)  - / p(t-t  ) h(t  ) dt 

™ -oo 

The  Fourier  Transform  of  above  equation  is 

P (f)  = P(f)  H(f ) 
m 

where  P(f),  H(f)  and  P (f)  are  the  Fourier  Transforms  of  p(t), 


h(t ) and  pni(t) , respectively.  The  magnitude  of  H(f)  can  be 
defined  as 


I H(f ) | - |Pm(f)|  / |P(f)| 


Applying  the  assumption  that  <P,,l!ff>  = 0,  we  can  derive  the  auto  spec- 
trum of  wave-induced  pressure  as  follows: 


* * * * * 
s *s  = <p  tt > • <rr  p > = <p  h rr> - <ir  ph> 

Pm^  npm  m m 


"M  A -k  ^ ^ 

= <P  TT>*  <n  P>  = <P  Tf>  -<TT  P>  = <P  P>  • <*T  TT > 


= S'- 
p 


(4.10) 


Therefore , 


S,(f)  “ S -S  / S„  « [ Co^  _(f)  + Qu^  _ (f )] / S (f ) 

Km  “m 


? Pmn  npm  ■ n • p n ' ' ' P„n  ' n 


(4.11) 


where  H (f)  H(f)  *=  |H(f)|^  ” 1,  which  was  shown  in  Figure  3.7. 


4.2.4  Cross  Spectrum  between  Wave  Height  and  Pressure  Fluctuation 
The  cross  spectra  between  wave  height  and  pressure  fluctua- 
tion can  be  u9ed  to  calculate  the  energy  and  momentum  transfer  from  wind 
to  waves  due  to  the  normal  pressure  acting  on  the  water  surface.  The 
cross  spectrum  between  wave  height  and  the  measured  pressure  signal  is 

S (f)  = <P*(f)rr(f)>  - Co  (f)  + i Qu  (f) 

Pmn  " Pm " Pm'1 


(4.12) 


while  the  cross  spectrum  between  wave  height  and  the  desired  pressure 
signal  is 


spn(f)  = <p*(f>n(f)>  = Spffin(f)  H(f)  = CoPn(f)  + 1 QuPn(f)  (4,13) 


where  H(f)  “ j H ( f ) | e » Cos  4> C f ) + Sin  Here  <f>(f)  is  the 

phase  angle  of  the  transfer  function  H(f)  and  was  determined  during 
calibration  of  the  pressure  measuring  system  as  shown  in  Figure  3.8. 


Substituting 


H(f  ) 


and  equation  (A. 12)  into  equation 


S (f)  * (Co  (f)  • Cos  4> ( f ) - Qu  (f)  • Sin  $(f)) 

pn  pn  pn 


+ i (Co  (f)  • Sin  <f>(f)  + Qu  (f)-Cos  <J>(f))  (A.1A' 

Pmn  p n 


By  comparing  the  real  and  imaginary  parts  of  equations  (A. 13)  and 
(A.1A),  the  cospectrum  and  quadrature  spectrum  between  the  wave  height 
and  pressure  fluctuation  are  seen  to  be 


Co  (f)  = Co  (f)  • Cos  4> ( f ) - Qu  (f ) • Sin  cf> ( f ) (A.l5) 

pn  p n p n 

m m 

Qu  (f)  = Co  (f)  • Sin  <fr(f)  + Qu  (f)-Cos  <j>(f)  (A. 16) 

pn  pmn  Pmn 

Based  on  equations  (A.  15)  and  (A. 16),  the  coherence  spectrum  and  phase 
between  wave  height  and  pressure  fluctuation  are 

RPn(f)  * [(Copn(f)  + Qupn(f))  ' <sp(f) *Sn(f ) ) ] 1/2  <«-17) 

©nn(f)  - tan'1  (Qu  ( f ) / Co  (f ) ) (A. 18) 

pn  pn  pn 
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Pm<f>  VO' 
H*(f ) * H(f ) i 


where  S (f)  - 
P 


<P  (f)P  (f)>  = S (f) 
m m Pm 


In  processing  any  functional  relationship  involoving  the  measured 
pressure  signal,  it  is  necessary  to  account  for  the  pressure  system 
response  as  indicated  by  the  results  of  calibration.  However,  system 
response  is  not  a factor  for  the  hot-film  data,  and  equations  (4.2), 
(4.3)  and  (4.4)  can  be  used  directly  for  calculating  cross  spectrum, 
coherence  spectrum  and  phase  difference  between  wave  height  and  velo- 
city fluctuations,  respectively. 


4.2.5  Energy  Transfer  from  Wind  to  Waves 

Since  it  is  the  pressure  components  out  of  phase  with  the 
progressive  waves  that  transfer  energy  and  momentum  from  the  wind  to 
the  wave  field,  the  quadrature  spectrum  between  wave  height  and  pres- 
sure in  equation  (4.16)  is  the  key  quantity  in  the  process.  By  sub- 
stituting equation  (4.16)  into  equations  (2.7)  and  (2.8),  the  energy 
and  momentum  transfer  from  wind  to  waves  due  to  normal  pressure  acting 
on  water  surface  can  be  expressed  as 


Sln(f)  - 2irf  *(Co  (f )•  Sin  #(f)  + Qu  (f).Cos  $(f))  (4.19) 

Pmn  V 

t (f)  - 4tt  2 f 2 • (Co  (f ) • Sin  <J.  (f ) + Qu  (f)-Cos  <J>  ( f ) ) / g (4.20) 

w 


It  can  be  seen  that  the  measured  pressure  fluctuation  p (t),  wave 

in 

height  n(t) , and  the  calibration  values  of  the  pressure  measuring 
H(f)  and  <Mf)  are  needed  for  such  a calculation. 


system 


4.2.6  Energy  Transfer  of  Nonlinear  Wave-Wave  Interactions 


Since  Barnett's  parameterized  equation  (2.17)  yields  satis- 
factory agreement  compared  with  Hasselmann's  theoretical  formulation, 
Barnett's  approximate  equations  are  used  in  this  study.  With  the  ex- 
perimentally determined  one-dimensional  energy  spectrum  E(f)  at  a 
certain  fetch  and  the  commonly  adopted  spreading  factor  d'(9)  for 
laboratory  data,  the  wave  parameters  e,  fQ  and  &o  can  be  calcula- 
ted. Consequently,  t and  ? can  be  evaluated  from  equations  (2.18) 

and  (2.19),  and  the  non-linear  wave-wave  interaction  S , (f)  can  be 

nl 

calculated  from  equation  (2.17).  Because  the  direction  of  wave  propa- 
gation is  from  shorter  fetch  to  longer  fetch,  the  information  for  waves 
at  shorter  fetch  could  be  considered  as  the  initial  condition  for  the 
rate  of  change  of  energy  between  two  consecutive  stations.  Therefore, 
the  wave  height  spectrum  E(f)  at  shorter  fetch  was  used  in  equations 
(2.17),  (2.21),  (4.22),  and  (4.23)  to  calculate  S (f),  Sd  g ( f ) , 

S ,(f),  and  S,  (f),  respectively, 
nl  ds 

It  is  to  be  noted  from  equation  (2.5)  that  the  net  source  function 
S(f)  can  be  calculated  from  the  experimentally  measured  wave  spectrum 
between  two  consecutive  stations  along  the  channel  and  can  be  expressed 

as 


S(f) 


V (f)  p g 

x w6 


AE(f) 

Ax 


32  g 

°-8 


w 


AE(f) 

Ax 


In  addition,  the  energy  input  to  the  wave  from  the  air  stream  S^tf) 
was  measured  and  can  be  calculated  from  equation  (4.19).  Consequently, 
the  experimental  results  of  the  sum  of  source  terms  S (f ) and 


S,  (f)  can  be  evaluated  as 
ds 


+ 5d.(f)  ■ S<f)  - Sln(f> 


(4.21) 


Unless  separate  and  independent  experimental  measurements  of  either 

Sni(f)  or  Sdg(f)  are  made,  it  is  not  possible  to  separate  these  two 

functions.  Futhermore,  it  seems  extremely  difficult  to  separate  these 

functions  experimentally.  In  order  to  calculate  the  non-linear  wave- 

wave  interaction  based  on  an  experimental  data,  S^s(f)  must  be 

assumed.  If  Hasselmann's  white-capping  dissipation  model  is  adopted 

? 

and  S.  (f)  “ -d(2irf)  E(f)p  g is  calculated  from  equation  (2.21), 
ds  w 

then  the  "experimental"  S . (f)  can  be  represented  by 

nl 

Sni<£)  - S(f)  - Sln(f>  - Sds(f)  (4.22) 


The  determination  of  the  constant  d 
the  next  section. 


in  S.  (f)  will  be  discussed  in 
ds 


4.2.7  Energy  Dissipation  of  White  Capping 

As  discussed  in  the  previous  section,  in  order  to  verify 
Hasselmann's  white  capping  dissipation  model,  the  source  term  Sn i ( f ) 
would  have  to  be  assumed  known.  If  Barnett's  approximate  representa- 
tion of  S , (f)  is  adopted,  the  calculated  or  "experimental"  white 
nl 

capping  dissipation  can  be  represented  as 


The  white  capping  dissipation  represented  by  equation  (2.21)  con- 
tains an  unknown  constant  d which  must  be  determined  from  the  experi- 
mental data.  In  the  high-frequency  equilibrium  range  of  a wave  spec- 
trum, the  net  source  term  S^(f)  must  be  zero.  Therefore,  we  can 
express  the  energy  and  mean  momentum  transfer  in  the  high  frequency 
range  as 


where  the  superscript  hf  indicates  the  high  frequency  regime. 

Since  the  energy  spectrum  in  the  equilibrium  range  can  be  approx- 
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imated  by  the  one-dimensional  form  ag  f , the  mean  momentum  transfer 
(in  the  x-direction)  due  to  white  capping  in  the  high-frequency  equili- 
brium range  can  then  be  represented  through  the  use  of  equation  (2.21) 
as  (Hasselmann,  1974) 


» (128  / 45tt)  d p g2a  16ir4  / 2*f  (4.25) 

w so 


■t 


TT  / 2 

where  1 / c-  (f)  = k / 2-rrf  = J (k  / 2irf)  F(f,6)  d0  / E(f) 

X X /«  x 


*/2 

= J (k  / 2itf)*Cos  6- (8/3ir)‘Cos  0*E(f)  de  / E(f) 

-til  2 


= (k  / 2irf ) (128  / 45tt)  = (2-irf  / g)  (128  / 45tt) 


Here  k^  Is  the  average  value  of  wave  number  k^  integrated  over  the 
normalized  directional  distribution  of  E(f),  is  the  water  den- 

sity, g is  the  gravitational  acceleration,  a is  Phillips  constant 
of  the  equilibrium  range  of  wave  spectrum,  and  f is  the  lowest 
frequency  at  which  the  net  energy  transfer  S(f)  equals  to  zero. 

Both  a and  f can  be  determined  by  the  experimental  results  of 

so 

wave  height  spectrum  and  net  energy  transfer  spectrum  S(f). 
Substituting  equation  (4.24)  into  equation  (4.25),  the  constant  d 
can  be  expressed  as 


d » 


45  n 
128 


2 n f 


so 


pwg  a 16w 


i -hf 
* (~  Tin 


-hf. 

\l> 


(4.26) 


__hf  x 

where  rV  “ f 2jrf  Qu  (f)  ttv  df  can  be  evaluated  from  the 
in  "r  pn  2irt 

r80 

experimental  value  of  Qu^  (f)  with  a chosen  f^. 

The  evaluation  of  the  term  needs  further  discussion.  The 

ni 

theoretical  results  of  non-linear  wave-wave  interaction  S , (f)  cal- 

ni 

culated  by  Barnett's  approximate  parametric  equation  seem  not  to  be 
adequate  in  the  high  frequency  range  [as  we  will  see  in  Figures (5 . 20) 
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through  (5.31)].  Some  modification  in  calculating  the  terms  can 

be  made,  based  upon  Hasselmann's  original  assumption  of  conservative 
wave-wave  interaction.  This  may  be  represented  as 


Anl(f)  « - J‘0S°Snl(f)  2if  df  + S Snl(f)  2?f  df  = 0 


(4.27) 


Consequently,  the  required  term  for  evaluating  d in  equation 

(4.26)  may  be  represented  as 


(4.28) 


After  d was  determined,  the  "theoretical"  results  for  S,  (f)  were 

ds 

calculated  by  equation  (2.21)  and  compared  with  "experimental"  results 
computed  by  equation  (4.23). 

4.2.8  Mean  Momentum  Transfer  from  Wind  to  Waves 

The  ratio  y defined  in  equation  (2.24)  represents  the 
magnitude  of  the  fraction  of  mean  momentum  transferred  to  waves 
divided  by  the  total  mean  momentum  transferred  across  the  wind-wave 
interface.  It  is  shown  in  Appendix  C that  the  mean  momentum  trans- 
ferred into  the  wave  field  can  be  expressed  as  [c.f.,  equation  ( C . 7 ) ] 
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and  the  total  mean  momentum  transferred  across  the  wind-wave  interface 
is,  [c.f.  equation  (C.6)] 

M + M = p + p u'^  - p u v 

w c 3x  a 3x  a 

In  order  to  estimate  the  value  of  y>  three  time-average  terms  in  the 

•>  o | f 

above  equation,  i.e.,  p 3n/3x,  u'  3q/3x  and  u v , have  to  be  eva- 
luated. The  data  p(t),  u'(t),  v’  (t)  and  n(t)  required  for  cal- 
culating these  three  terms  are  in  the  time  domain.  Since  it  is  ex- 
tremely difficult  to  deduce  the  wave-induced  (tilda-term)  and  background 
turbulent  (prime-term)  terms  from  the  total  signal  for  the  wind  wave 
spectrum,  an  alternative  way  was  employed  to  calculate  these  three 
terms.  Dobson  (1971)  showed  that  the  time-average  results  of  these 
terms  equals  the  real  part  of  the  integration  of  their  cross  spectra 
over  the  appropriate  frequency  range.  Therefore,  the  cross  spectra 
between  p and  3n/3x  (or  p and  n)»  u'^  and  3n/3x  (or  u'^  and 
q),  and  u'  and  v'  (or  u and  v)  should  be  calculated  first  and  then 
integrated  between  certain  lower  and  upper  frequency  limits.  The  lower 
limit  of  the  integration  is  0.078  Hz  which  is  the  bandwidth  of  the 
corresponding  spectral  estimates  calculated  from  the  relation  Af  = 

1/ (rN)  • 1/(0.025  x 512)  - 0.078  Hz,  where  T (0.025  sec)  is  the  sam- 
pling rate  and  N(512)  is  the  number  of  data  for  each  block  of  the 
Fourier  Transform.  The  upper  limit  is  10.0  Hz  which  is  the  upper  limit 
for  a gravity  wind-wave.  The  cross  spectra  associated  with  the  deriva- 
tive term  3q/3x  can  be  simplified  to  that  of  n only  by  multiplying 
a factor  of  (-  ik)  as  a consequence  of  the  form  of  the  Fourier 


Transform  of  3q/3x.  Because  of  appearence  of  i («  /-T) , the  quadra- 
ture spectrum  iQu(f)  was  adopted  in  order  to  obtain  the  real-value 
results. 


The  detailed  expressions  for  calculating  the  three  key  terms  were 
then  formulated  as  follows.  The  first  term  p 3n/3x  can  be  expressed 


p U - S (-  ik>  ti  Qupn(f)]  df 


E ki  Qupn(fi)  Af 

f "0.078  V 


where  kj  = (2Trf-i)  / g is  the  wave  number,  Qu  _,(f)  is  the  quadra- 

P '1 

ture  spectrum  between  pressure  and  wave  height,  and  Af  = 0.078  Hz  is 
the  bandwith  of  the  spectrum  estimates.  The  more  explicit  form  of 
p 3n/ 3x  can  also  be  written  as 


-*=■  10-°  7 2 

E (*»  fj.  / g)  Qunn(fi>  Af 

3X  ft-0.078  1 Pn 


The  quadrature  spectrum  Qu  (f^)  can  be  calculated  based  on  equation 
(A. 16). 

The  calculation  of  the  second  term,  u'2  3n/3x  involved  extreme 
difficulty  in  calculating  the  Fourier  component  of  u'^,  i.e.  U'U'(f), 
so  the  following  relation  was  applied  to  simplify  the  calculation: 


(u  - u)2  - u2  * (u  + u ' ) 2 - u2  * u2  + 2uu'  + u'2  - u2  = u'2 


i v l .** 


where  the  total  signal  u equal  (u  + u + u’)  and  the  quantities  with 
the  bar,  tllda  and  prime  represent  mean,  wave-induced,  and  background 
turbulent  quantities  respectively.  Through  the  use  of  the  above  rela- 
tions, the  second  term  (In  the  time  domain)  can  be  rewritten  as 


( — .2  3q  - 2 9n 

p (u  - u)  — - p u — 
a 9x  a 9x 


Based  on  the  previous  discussion,  these  two  terms  can  be  expressed  in 
the  frequency  domain  in  a form  similar  to  equation  (4.29),  viz.. 


— 2 10.0  « o 

p.<u  - 10  S?  - pa  fi.J  078  <4’  fl  ' 8) 


2 Sn  iU>U  0 0 

Pa  5 3?  “ pa  E (4ir  fl  ' g)  Qu  9 

3 8x  3 fi-0.078  1 u2n  1 


• j a 10 .0  2 9 

Pa  u'  9?  = pa  fiJ>078  (4*  fi  ' 8)  Af[Qu(u-a)2n(fl) 


(4.30) 


2 

where  anc^  Qu~2  are  Quadrature  spectra  of  (u  - u)  and 

-2  U ^ 

u with  respect  to  n,  respectively.  Since  the  time  average  quantity  IT 
can  be  easily  obtained,  Qu^  xr)^n  can  ca^culated  from  equation 
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(4.2)  by  replacing  x and  y by  (u  - u)  and  n respectively. 

Because  of  the  difficulties  involved  in  decomposing  the  wave  induced 

signal  u(t)  (or  U(f)),  related  to  the  wave  spectrum,  from  the  total 

signal  u(t),  the  quantity  Qu,„  can  not  be  calculated  directly  from 

u^n 

equation  (4.2).  Therefore,  the  assumption  described  by  equation  (4.9) 

in  Section  4.2.2,  i.e.,  U(f)  = <U(f)  TT*(f)>  TT(f ) / <TT*(f)  1T(f)>  is 

-2 

employed  to  calculate  the  Fourier  Transform  of  u (t) , denoted  as  UU(f), 

~2  * 

then  the  cross  spectrum  between  u (t)  and  n(t),  i.e.,  [UU(f)]  TT(f), 

and  its  quadrature  spectrum  Qu  „ (f).  The  derivation  are  shown  below. 

Q^n 

If  a(t)  and  b(t)  are  any  two  time  signals  and  A(f)  and  B(f) 
are  their  respective  Fourier  Transforms,  then  the  Fourier  Transform  of 
the  convolution  of  a(t)  and  b(t)  is  (Bracewell,  1965) 

F{a (t)  * b(t)}  = A(f ) B(f ) 


The  convolution  of  the  Fourier  Transform  of  a(t)  and  b(t)  is 

A(f ) * B(f ) = F(a(t)  b (t) } 

where  * denotes  convolution  and  F { } represents  Fourier  Transform. 
The  definition  of  convolution  is 


or 


a(t)  * b(t)  - J a(t-t1)  b(t1)  d^ 


A(f)  * B(f ) - J A(f-f1)  B(fx)  dfx 
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Based  on  the  definition  of  convolution  and  the  relationship  described 

2 

above,  we  can  express  the  Fourier  transform  of  u (t)  as 


UU(f ) ■=  F{u(t)  u(t) } = U(f)  * U(f)  = J U(f-f1)  U(fx)  dfx 


or 


10.0 

UU(f  ) = T.  U(f  -f  ) U(f.)  Af 
f =-10.0  1 2 2 


( 4.31) 


where  i = 1,  128,  f = 0.078  Hz,  f.,,  = 10.0  Hz,  Af  = 

1 1^0 


0.078  Hz 


The  procedure  for  calculating  Qu_2  (f)  can  be  summarized  as 

u^n 

follows.  The  spectral  average  <U(f)  TT(f)>,  Fourier  components  iT(f) 

and  the  average  wave  spectrum  <TT*(f)  TT(f)>  are  calculated  first  in 

order  to  approximate  the  Fourier  components  of  wave-induced  velocity 
**  - ? 

U(f).  Then  the  Fourier  component  of  u (t),  i.e.,  UU(f),  can  be  eva- 
luated  by  equation  (4.31)  and  hence  the  cross  spectrum  between  u^(t) 

a,  « ^ 

and  ri(t),  i.e.,  UU  (f)  fl(f).  Finally,  the  quadrature  spectrum 
Qu-9  (f)  is  determined  by  the  imaginary  part  of  the  cross  spectrum 

.J{  " 

UU  (f)  TT( f ) . After  determining  Qu,  _.2  (f ) and  Qu  „ (f),  the  mean 

' u u ) 0 u^-p 

momentum  transfer  from  wind  to  waves  due  to  the  wave-induced  turbulent 


stresses  p u’2  3n/9x  can  be  calculated  by  equation  (4.30). 

a 

/ 

The  third  term  -p  uTvr  can  be  rewritten  by  applying  the  pro- 
perties of  u v'  * 0 and  v u'  = 0 and  the  relations  u = IT  + u + u' 
and  v - v + v + v',  namely, 


p u'v’ 

a 


p^Ku  + u')  (v  + v')  - uv] 


- p [ (u  - u)  (v  - v)  - uv] 


(4.32) 
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The  total  Reynolds  stresses  (u  - U)  (v  - Y7)  can  be  calculated 

by  multiplying  the  time  signals  of  u(t)  and  v(t)  after  removing 

their  mean  values.  The  wave-induced  Reynolds  stresses  uv  can  be 

evaluated  by  integrating  the  real  part  of  the  cross  spectrum  of  u and 
★ 

v,  i.e.,  Re{U  (f)  V ( f ) } or  Co-^(f),  over  the  interested  frequency 
range.  Thus, 


uv  = J Re{ U (f)  V(f ) } df 


10.0 

l 

f 1=0.078 


Co~ - (f j ) &f 

UV  1 


(4.33) 


where  Co~~(f)  is  the  cospectrum  of  U (f)  V(f).  Similarly,  the 
assumptions  of  U(f)  « <U(f)  fl  (f)>  ff(f)  / <TT*(f)  TT(f)>  and 
V(f)  = <V(f)  IT  (f)>  TT ( f ) / <TT  (f)  TI(f)>  are  employed  to  calculate 
U*(f)  V(f). 

The  mean  momentum  transferred  into  waves  M can  be  determined 

w 

by  the  sum  of  equations  (4.29)  and  (4.30).  The  total  mean  momentum 

transferred  across  the  wind-wave  interface  M + M can  be  determined 

w c 

by  the  sum  of  equations  (4.29),  (4.30)  and  (4.32).  Finally,  the  value 

Y can  be  evaluated  by  taking  the  ratio  of  M over  (M  + M ). 

w w c 

4.2.9  Mean  Velocity  Profile  and  Friction  Velocity 

The  mean  velocity  profiles  were  calculated  by  equation  (3.1) 
using  the  averaged  outputs  of  the  pitot-static  tube.  Many  experimental 
investigators  in  both  laboratory  and  field  have  shown  that  the  mean 


-49- 


■•XT 

• -O*'  * 


velocity  profile  u/u*  is  logarithmic  with  u^y  / v and  may  be  repre- 
sented by 


u* 


1 


Jin 


+ c 


(4.34) 


where  y is  the  distance  measured  from  mean  water  surface,  k = 0.4 
is  von  Karman  constant,  and  v is  the  kinematic  viscosity  of  air.  The 
friction  velocity  and  constant  c were  determined  by  a least 

square  fit  from  the  experimental  data.  The  measured  velocity  profiles 
are  shown  in  Figures  4.1  and  4.2.  The  velocity  profiles  near  the  water 
surface  deviate  from  the  assumed  logarithmic  distribution,  due  mainly 
to  the  effects  of  the  waves.  Similar  results  were  also  found  by  Wu 
(1968),  Sutherland  (1967)  and  Chao,  et  al.  (1976).  Consequently,  the 
data  in  this  region  were  excluded  in  the  least  square  fit.  Only  the 
data  in  the  central  portion  of  the  profiles  were  used  to  determine  u*. 


4.2.10  Equilibrium  Constant,  Spectral  Peak  Frequency,  and  Total 
Energy  of  Wave  Spectrum 

Phillips  (1958)  found  that  the  functional  form  of  the  wave 
spectrum  in  the  equilibrium  range  can  be  represented  as 

E(f)  ■ a g^  f or  E(tj)  = a g^  w (4.35) 


i 4 

where  a - a(2ir)  is  a constant,  g is  gravitational  acceleration, 
and  u “ 2n  f . 

The  constant  a in  this  study  was  calculated  by  two  different 
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methods.  The  first  method  is  to  use  the  least  square  fit  of  the 


measured  wave  height  spectrum  E(f)  to  equation  (4.35)  in  a given 

frequency  range.  Pierson  (1959)  and  Sutherland  (1967)  found  that  the 

appropriate  frequency  range  is  that  lying  between  the  second  harmonic 

of  the  dominant  wave  frequency  and  the  highest  limiting  frequency  of 

the  gravity  wave  spectrum.  The  maximum  frequency  of  a gravity  wave 

spectrum  in  this  study  is  about  10  Hz.  The  second  method  is  based  on 

2 -5 

an  average  technique.  Since  a equals  E(f)  / (g  f ),  a series  of 
a values  can  be  obtained  as  a function  of  frequency  f in  the  equi- 
librium range  of  a given  gravity  wave  spectrum.  Then  an  average  value 
of  a can  be  obtained  over  the  range  of  the  frequencies.  The  results 
obtained  from  both  techniques  are  presented  in  next  chapter. 

The  spectral  peak  frequency  of  each  wave  spectrum  is  the  frequency 
at  which  the  spectral  density  is  a maximum.  A simple  relationship  that 
exists  between  the  peak  frequency  f and  fetch  x may  be  expressed 
as 


f 

m 


(4.36) 


where  f * f u.  / g is  the  nondimensional  peak  frequency  and 
ro  m * 

x » x g / u^  is  the  nondimensional  fetch.  The  constant  c^  and  the 

exponent  n were  determined  by  a least  squares  fit  of  the  experimental 
data. 

Similarly,  a simple  relationship  between  the  total  energy 
00 

e = J E(f)  df  in  a wave  spectrum  and  the  fetch  also  exists  and  can  be 
o 

represented  in  nondimensional  form  as 
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-m 

E = X 

A 2 4 

where  e * e g / is  the  nondimensional  total  energy  of  the  wave 
spectrum  and  x is  the  nondimensional  fetch.  Again  the  constant  c0 
and  the  exponent  m were  determined  from  the  experimental  data  by  a 
least  squares  fit. 


5.  PRESENTATION,  COMPARISON,  AND  DISCUSSION  OF  RESULTS 


In  this  chapter,  the  general  characteristics  of  the  wave,  wind 
and  wave-induced  fields  are  presented  first.  Then,  the  major  subjects 
of  this  study  including  the  energy  and  momentum  transfer  from  wind  to 
waves,  nonlinear  wave-wave  interactions  and  white-capping  dissipation 
are  given.  The  comparison  and  discussion  of  the  results  are  also 
provided  in  each  section. 

5. 1 Characteristics  of  the  Wave  Field 
5.1.1  Wave  Height  Spectrum 

The  measured  power  spectra  of  wave  height  at  five  stations 

along  the  channel  are  shown  in  Figure  5.1  (a),  (b) , (c)  for  three  wind 

velocities  7.09,  8.01  and  8.88  m/sec.  The  mean  wind  velocities  were 

measured  at  a fetch  of  12.61  m.  All  the  spectra  have  the  common 

characteristics  of  a steep  forward  face,  a sharp  peak  at  the  dominant 

frequency,  and  the  second  hump  at  the  second  harmonic. 

Since  the  waves  are  growing  along  the  channel,  stations  at  larger 

fetch  have  a higher  spectral  peak  implying  more  energy  per  frequency 

band.  The  normalized  total  wave  energy  e for  three  wind  speeds  as 

a function  of  the  normalized  fetch  x is  presented  in  Figure  5.2  (a). 

The  least  square  fitted  curve  indicates  the  approximately  linear 

relationship  between  the  normalized  total  wave  energy  and  the  nor- 

/>  .1.23 

malized  fetch,  i.e.,  e 11  x . This  curve  is  comparable  witii  other 
measurements  carried  out  both  in  the  ocean  and  laboratory.  Our  experi- 
mental results  were  also  plotted  in  Figure  5.3  (b)  whose  form  was 
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introduced  by  Hasselmann,  et  al.  (1973),  for  comparison. 

The  spectral  peak  frequencies  shift  from  high  to  low  frequency 
along  the  channel  as  shown  in  Figure  5.1.  Similar  behavior  also  occurs 
in  the  ocean  wave  spectrum  measured  at  a fixed  location  with  increasing 
wind  or  measured  along  the  fetch  in  the  direction  of  wave  propagation 
for  a fixed  wind  speed.  This  effect  is  shown  later  to  be  caused  by  the 
nonlinear  wave-wave  interaction  mechanism  which  transfers  energy  from 
the  intermediate  frequency  to  the  lower  frequency  range.  The  non- 
dimensional  spectral  peaks  f as  a function  of  nondimensional  fetch 
x are  presented  in  Figure  5.2  (b) . Figure  5.3  (c)  shows  also  the 
comparison  of  our  laboratory  data  for  f^  vs.  x with  field  data. 

Most  of  the  results  indicate  that  the  normalized  spectral  peak  f 

is  proportional  to  the  -0.33  power  of  the  normalized  fetch  x,  i.e., 

.-0.33 

f “ x . The  present  data  can  be  also  fitted  to  this  relationship 

A “0.44 

very  well,  although  the  best  curve  fit  results  in  f x 

Finally,  according  to  the  previous  and  present  results  of  e,  f and 

x,  the  fetch  dependence  of  f and  e is  established. 
r m 

The  variations  of  some  spectral  components  (1.5  Hz  to  10  Hz)  of 
the  wave  at  five  stations  are  shown  in  Figure  5.4  for  wind  speeds  of 
7.09,  8.01  and  8.88  m/sec.  The  graphs  indicate  that  the  frequency 
components  equal  to  or  less  than  2.5  Hz  are  still  growing  and  not  yet 
saturated.  The  components  larger  than  2.5  Hz  are  sometimes  growing 
and  sometimes  decaying,  which  reveals  the  exchange  of  energy  among 
different  wave  components  associated  with  the  nonlinear  wave-wave 
interactions.  Not  much  variation  was  found  for  the  frequency  compo- 
nents larger  than  10  Hz  which  is  in  the  capillary  wave  range. 
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5.1.2  Equilibrium  Constant  of  Wave  Height  Spectrum 


The  growth  of  waves  under  the  action  of  wind  can  not 
continue  indefinitely.  The  existing  experimental  data  for  wind  wave 
spectra  obtained  in  both  the  laboratory  and  the  field  indicate  that 
there  exists  a "saturated"  or  equilibrium  range  at  high  frequency. 
Phillips  (1966)  argued  that  there  are  two  possible  physical  reasons 
which  prohibit  the  growth  of  the  wave  spectrum  and  lead  it  to  an 
equilibrium  state.  First,  if  the  waves  become  too  steep,  local 
instability  causes  the  wave  to  break.  Consequently,  energy  is  lost 
from  the  wave  field.  The  appearance  of  intermittent  patches  of  white 
caps  on  the  ocean  surface  indicates  the  local  instability  leading  to 
wave  breaking.  Another  mechanism  for  the  limitation  of  wave  growth 
arises  from  the  formation  of  capillary  waves  in  the  front  of  steep 
primary  waves.  The  lower  limit  of  the  frequency  of  a capillary  wave 
is  about  10  Hz.  Therefore,  the  frequency  range  of  the  equilibrium 
portion  of  a wave  spectrum  is  between  the  frequency  of  the  second 
harmonic  in  the  wave  spectrum  (Pierson,  1959)  and  10  Hz.  Figure  5.1 
shows  the  wave  spectrum  at  five  stations  along  the  channel  for  three 
wind  speeds  of  7.09,  8.01  and  8.88  m/sec  and  indicates  the  existence 
of  this  equilibrium  range. 

Based  upon  purely  dimensional  consideration,  Phillips  (1958) 
proposed  equation  (4.35)  representing  the  one-dimensional  equilibrium 
range  spectrum.  Furthermore,  Phillips  found  a to  be  a universal 
constant  based  on  the  existing  wave  spectra  data.  On  the  other  hand, 
recent  JONSWAP  ocean  spectrum  data  (Hasselmann,  et  al , 1973)  indi- 
cated that  the  equilibrium  range  of  an  ocean  spectrum  is  represented 
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by  the  energy  balance  between  the  energy  Input  from  turbulent  air, 
the  white  capping  or  wave  breaking,  and  nonlinear  wave-wave  inter- 
actions. Hasselmann  (1974)  proposed  the  dissipation  by  white  capping 
1 1 

as  S (k)  = - d (2irf)  E(k)  , and  found  the  equilibrium  range  of  a 
wave  spectrum  can  be  approximated  by  f~5. 

The  values  of  a in  the  equilibrium  range  wave  spectrum  were 
determined  by  two  different  methods  as  discussed  in  Section  4.2.10. 
Figure  5.2  (c)  shows  the  variation  of  a as  a function  of  non- 
dimensional  fetch,  x . It  indicates  that  a is  not  a universal 

constant  and  is  definitely  fetch  dependent.  The  experimental  results 

.-0.65 

indicate  a “ x . The  same  set  of  data  (a  vs.  x)  was  replotted 
in  Figure  5.3  (a)  for  comparison  with  the  results  of  other  investi- 
gators. This  set  of  data  is  in  the  region  of  small  non-dimensional 
fetch  and  compares  quite  favorably  with  other  existing  data. 

From  the  results  of  Section  5.1.1  and  this  section,  we  find  the 
nondimensional  peak  frequency  £ , nondimensional  total  wave  energy 

e , and  Phillips  equilibrium  constant  a to  be  fetch  dependent. 

The  best  curve  fit  for  these  variables  with  respect  to  the  nondi- 
mensional fetch  y yield 

- . -0.44  ' - 1.23  . -0.65 

f a x , e “ x and  a x 

m 

~ * — 4 

For  a self-similar  spectrum  e ~ af^  ; accordingly,  the  rela- 
tionship among  f^  , e and  a can  be  constructed  and  the  exponents 

nf  , n and  n in  the  power  laws  for  f , e and  a should 
tea  r m 

satisfy  the  following  relation: 
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4n,.  - n + n =0 

f a e 

Here  n,  - -0.44  , n * -0.65  , and  n = 1.23  for  the  present 
f a £ r 

experimental  results.  Substituting  these  values  into  the  above  rela- 
tion yields 


4 (-0.44)  - (-0.65)  + 1.23  = 0.12 

The  resultant  value  0.12  lies  within  the  error  of  about  ±0.24  esti- 
mated from  the  individual  error  band  of  each  curve  fitted  exponent. 
Finally,  based  on  the  JONSWAP's  (Hasselmann,  et  al.  1973)  and  the 
present  results  of  the  relationships  of  f , e and  a to  x , the 
fetch  dependence  of  ? , e and  a is  again  established. 

5 . 2 Characteristics  of  the  Wind  Field 

The  mean  velocity  profile  and  the  power  spectra  of  the  fluc- 
tuating velocity  components  u'  and  v'  characterize  the  mean  and 
fluctuating  wind  fields.  The  mean  air  velocity  profiles  were 
presented  in  Figures  4.1  and  4.2.  The  vertical  line  in  the  upper 
part  of  the  profile  is  the  free  stream  region.  The  central  part 
represents  the  logarithmic  region  of  a typical  turbulent  velocity 
profile.  The  lower  region  of  the  velocity  profile  close  to  the  air- 
water  interface  deviates  from  the  logarithmic  profile.  There  the 
velocity  magnitudes  are  larger  than  those  represented  by  a logarith- 
mic law.  The  deviation  from  a conventional  logarithmic  profile 
and  the  increase  in  magnitude  may  be  caused  by  the  generation  of  the 
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surface  drift  current.  Consequently,  the  friction  velocity  u*  was 
evaluated  from  the  best  fit  of  equation  (4.34)  based  on  the  experi- 
mental data  in  the  logarithmic  region.  The  calculated  results  of 
u*  are  tabulated  in  Table  (5.1). 

Hie  power  spectra  of  the  turbulent  velocity  components  u'  and 
v'  were  evaluated  from  the  experimental  data  taken  closest  to  the 
air-water  interface  (i.e.,  0.508  cm  above  the  highest  wave)  and 
sampled  at  a rate  of  5000  Hz  for  a duration  of  five  minutes.  The 
spectra  were  normalized  with  respect  to  their  respective  variance 
and  are  shown  in  Figures  5.5  to  5.8.  Figures  5.5  (a),  (b),  and  (c) 
and  5.6  (a),  (b) , and  (c)  represent  the  u'  and  v'  normalized 
power  spectra  for  three  wind  speeds  at  the  two  fetches  of  6.51  m and 
15.66  m,  respectively.  The  u'-component  is  consistently  larger 
than  the  v' -component  in  the  lower  frequency  range.  However,  for 
frequencies  higher  than  approximately  20  Hz,  the  trend  is  reversed. 
Figures  5.7  (a)  and  (b)  show  the  comparison  of  the  normalized  u' 
and  v'  spectra  for  three  different  wind  speeds  at  the  two  fixed 
fetches  of  6.51  m and  15.66  m,  respectively.  Both  the  u'  and  v' 
spectra  are  independent  of  wind  speed  in  the  inertia  subrange  where 
the  -5/3-law  holds.  However,  such  is  not  the  case  in  the  high 
frequency  range.  In  order  to  reveal  more  details  of  the  spectra 
variation  in  the  high  frequency  range,  the  u'  and  v'  component 
spectra  were  replotted  for  a fixed  wind  speed  at  5 different  fetches 
as  shown  in  Figures  5.8  (a),  (b),  and  (c).  The  spectra  indicate 
that  there  is  a consistent  fetch  dependence  for  all  wind  speeds 


According  to  a discussion  of  the  general  behavior  of  a turbu- 
lence power  spectrum  by  Hlnze  (1975),  the  low  frequency  portion  of 
a turbulence  spectrum  is  governed  by  the  characteristics  of  the  large 
eddies  and  is  dependent  on  the  conditions  of  the  formation  of  the 
spectrum.  On  the  other  hand,  the  high  frequency  portion  of  a 
turbulence  power  spectrum  is  governed  by  the  turbulence  energy  and 
viscous  dissipation  and  is  independent  of  the  condition  of  the  for- 
mation of  the  spectrum.  The  normalized  power  spectra  presented  in 
Figures  5.5  to  5.7  are  used  to  study  the  direct  wave  effect  on  the 
turbulence.  As  indicated  in  Figures  5.5  to  5.7,  the  u' -component 
spectrum  seems  not  to  be  affected  by  the  wave.  The  spectrum  shape 
of  u'  is  proportional  to  f-1  and  is  consistent  with  the  results 
of  Tchen's  (1953)  study  of  turbulent  kinetic  energy  in  a shear  flow. 
The  eddy  diffusivity  determined  by  larger  eddies  might  be  regarded 
as  a governing  parameter  in  determining  the  characteristics  of  the 
turbulence  spectrum.  Upon  examination  of  the  wind  wave  spectrum  as 
shown  in  Figure  5.1,  the  range  of  the  wave  frequency  (f  < 10  Hz) 
is  in  the  low  frequency  portion  of  the  turbulent  spectrum.  The 
non-dimensional  v' -component  spectrum  deviates  consistently  from 
the  f“ 1 law,  and  a spectral  peak  appears  around  the  dominant  wave 
frequency.  It  shows  the  strong  wave  effect  on  the  v' -component 
spectrum  and  reveals  the  significant  contribution  of  the  wave  to  the 
formation  of  v'-component  spectrum  in  the  low  frequency  range. 
Similar  results  were  obtained  by  Young,  et  al.  (1973)  and  McIntosh, 
et  al.  (1975). 

In  e inertial  subrange  where  the  spectra  obey  the  f-5/3  law, 


the  total  energy  supply  to  the  turbulence  is  practically  equal  to 
the  total  viscous  dissipation.  It  is  clear  in  Figures  5.7  (a)  and 
(b)  that  such  equilibrium  conditions  were  established  for  both  u' 
and  v'  components  (20  Hz  < f < 500  Hz)  along  the  channel  [6.51  m, 
short  fetch;  15.66  m,  largest  fetch]  at  all  three  wind  speeds.  In 
addition,  the  existence  of  the  log-region  of  the  mean  velocity 
profile  as  shown  in  Figures  4.1  and  4.2  is  consistent  with  the 
existence  of  the  inertial  subrange  of  the  u'  and  v'  turbulence 
power  spectra.  For  frequencies  higher  than  those  of  inertial  sub- 
range, the  power  spectra  decays  as  f~ 1 if  the  viscous  dissipation 
is  the  only  governing  parameter.  It  can  be  seen  clearly  from 
Figures  5.8  (a),  (b) , and  (c)  that  the  measured  spectra  have  not 
reached  the  viscous  dissipation  range  as  evidenced  by  their  spectra 
decay  rate.  The  measured  spectra  depend  on  fetch  at  all  wind  speeds. 

5 . 3 Characteristics  of  the  Wave-Induced  Field 

The  wave-induced  field  is  represented  by  the  power  spectra  of 
the  wave-induced  velocity  components  u and  v and  the  wave-induced 
pressure  p , as  well  as  the  coherence  and  phase  difference  between 
wave  height  n and  these  wave-induced  quantities.  The  power  spectra 
of  wave-induced  velocity  components  u and  v and  the  wav-induced 
pressure  are  calculated  from  equations  (4.8)  and  (4.11)  respectively. 
The  coherence  and  phase  relations  between  n and  u and  v are 
obtained  from  equations  (4.3)  and  (4.4),  whereas  the  coherence  and 
phase  relations  between  n and  p are  calculated  by  equations 
(4.17)  and  (4.18).  Because  all  of  the  plots  are  similar,  typical 


sets  of  these  data  taken  at  0.508  cm  above  the  highest  waves  at  a 
short  fetch  of  6.51  m and  at  the  longest  fetch  of  15.66  m for  three 
wind  speeds  of  7.09,  8.01  and  8.88  m/sec  are  shown  in  Figures  5.9  to 

5.14. 

The  power  spectra  of  wave-induced  pressure  and  fluctuating 
velocities  have  shapes  similar  to  that  of  the  wave  height  spectrum 
as  shown  in  parts  (a)  and  (b)  of  these  figures.  Therefore,  the  flow- 
field  near  the  wave  surface  is  strongly  coupled  to  the  interfacial 
wave  motion. 

The  first  and  second  peaks  in  the  coherence  spectra,  as  shown 
in  parts  (c),  (d) , and  (e)  of  these  figures,  are  associated  with  the 
corresponding  harmonics  of  the  wave  spectrum.  These  indicate  the 
strong  coupling  between  the  wave-induced  field  and  wave  motion  at  the 
dominant  wave  frequency  (first  harmonic)  and  its  second  harmonic. 

In  addition,  the  high  values  of  coherence  imply  that  the  measured 
phase  relationships  between  the  wave-induced  quantities  and  the  wave 
motion  are  reliable.  As  shown  in  parts  (c) , (d)  and  (e)  of  these 
figures,  the  phase  angles  are  nearly  constant  near  the  wave  dominant 
frequencies.  The  phase  angles  scatter  beyond  the  dominant  fre- 
quencies. The  weak  interactions  between  the  two  flowfields  may  be 
responsible  for  the  phase  angle  scatter.  The  coherences  of  p,  u 
and  v with  respect  to  q ranged  from  0.70  to  0.95,  0.54  to  0.84 
and  0.87  to  0.95,  respectively,  at  the  dominant  frequency  along  the 
channel  for  all  cases.  The  plots  of  coherences  R against  non- 
dimensional  fetch  x are  shown  in  parts  (a),  (d)  and  (e)  of  Figure 

5.15.  It  was  found  that  the  coherences  between  n and  p and  q 
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and  u Increase  very  slightly  with  increasing  x,  whereas  the  co- 
herences between  n and  v are  almost  fetcli  independent.  The 
coherences  of  v and  p with  respect  to  n are  larger  than  that 
of  u with  respect  to  p . It  can  hardly  be  explained  by  the 
present  data.  One  possible  reason  may  be  due  to  the  effect  of  the 
interactions  among  the  mean  air,  wave-induced  and  turbulence  fields. 
Further  study  about  these  interaction  fields  is  needed  in  order  to 
fully  understand  the  dynamic  structure  of  the  air-sea  interface. 

It  can  be  seen  from  the  parts  (c),  (d) , and  (e)  of  Figures 
5.9  to  5.14  that  the  phase  angles  between  p,  u and  v with 
respect  to  n are  consistent  along  the  channel  (from  x = 6.51  m 


to  X = 

15. 

66  m) 

at 

all 

wind  speeds,  namely 

p leads  p 

by  220° 

to  230° 

(or 

P 

lags 

n 

by  130°  to  140°) , 

u leads  p 

by  10°  to 

o 

o 

/—V 

o 

u 

lags 

n 

by 

320°  to  350°)  and 

p leads  v 

by  65°  to 

75°  (or 

n 

lags 

V 

by 

285°  to  295°).  The 

experimental 

data  for 

the  phase  angles  between  p and  p,  u and  v at  the  dominant 
wave  frequency  along  the  channel  were  shown  in  parts  (b) , (c)  and  (f) 
of  Figure  5.15.  It  seems  that  the  phase  angles  between  wave  and 
wave-induced  fields  are  independent  of  the  nondimensional  fetch,  x. 
According  to  potential  flow  theory,  the  fluctuating  pressure  at  the 
water  surface  of  sinusodial  waves  propagating  in  the  same  direction 
as  a uniform  wind  field  is  (Lamb  1932,  §232) 


P“_Pa8n  [1  + (1  - u,  / c)2  ] 


where  is  the  mean  equivalent  uniform  wind  speed  and  c is  the 
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phase  speed  of  the  waves.  This  expression  predicts  a 180°  phase 
difference  between  n and  p for  a potential  flow.  An  experi- 
mental data  for  the  phase  angle  between  n and  p agrees  well 
with  the  field  data  obtained  by  Dobson  (1971)  and  Elliott  (1972)  , 
i.e.,  p leads  n by  225°  on  the  average.  Figure  5.16  shows 
the  comparison  of  the  phase  angle  data  of  p with  existing  experi- 
mental results  obtained  both  in  the  field  and  laboratory.  Although 
the  data  are  scattered,  a data  trend  was  found.  It  indicates  that 
the  phase  differences  between  r\  and  p approach  180°  for  large 
c/u*  (lower  wind  speed)  and  tend  to  deviate  from  180°  for  smaller 
c/uA  (higher  wind  speed).  The  difficulty  in  obtaining  accurate 
experimental  phase  angles  and  the  possible  effect  of  variable  tur- 
bulence characteristics  on  wave-induced  field  may  be  responsible 
for  the  scatter  of  the  experimental  data  in  the  higher  wind  speed 
range.  The  phase  angle  between  n and  u and  n and  v close 
to  wave  surface  could  be  derived  from  the  boundary  conditions  at 
the  wave  surface  (Yu,  et  al.,1973).  The  results  are  that  n leads 
u and  v by  180°  and  270°,  respectively,  inside  the  critical  layer 
(i.e.  , small  values  of  u^/c)  where  the  mean  wind  velocity  equals 
the  phase  speed  of  wave,  whereas  q leads  u and  v by  0°  and 
90°,  respectively,  outside  the  critical  layer  (i.e.,  larger  values 
of  Uro/c)  (Yu,  et  al. , 1973).  Since  the  wind  speeds  of  our  mea- 
surement are  relatively  high,  all  the  experimental  results  are  in  the 
range  of  large  u<r,/c.  The  present  phase  angle  results  are  plotted 
in  Figure  5.17  and  compared  with  those  of  Yu,  et  al. , 1973.  Although 
Yu's  data  was  obtained  in  a wave-following  system  over 
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artificially-generated  simple  progressive  water  waves,  there  is  good 
agreement  seen  if  Yu's  data  are  extended  to  the  region  of  large 
Uod/c.  Again,  the  results  deviate  from  the  theoretical  predictions 
based  on  the  boundary  condition  at  the  air-water  interface  (i.e., 
n leads  u and  v by  0°  and  90°  for  large  Uoo/c  ratio). 
Experimental  difficulties  preclude  us  from  measuring  true  pressures 
at  the  air-water  interface.  The  distance  of  the  pressure  probe  to 
the  interface  may  cause  such  deviations  from  the  theoretical  esti- 
mates. 

5 . 4 Energy  Transfer  from  Wind  to  Waves 

The  energy  transfer  from  wind  to  waves  is  through  the  work 
done  by  the  normal  wave-induced  pressure  p,  wave  induced  turbulent 
Reynolds  stresses  u'v'  and  v'v',  and  the  mean  turbulent  Reynolds 
stresses  u'v’  acting  on  wave  surface  as  shown  in  equation  (D.12) 
of  Appendix  D.  The  results  of  calculation  based  on  the  experimental 
data  were  tabulated  in  Table  (D.l).  The  energy  transferred  to  the 
wave  field  due  to  the  normal  pressure  term  is  calculated  from  equa- 
tion (4.1*)  and  all  other  terms  are  calculated  by  employing  a tech- 
nique similar  to  that  for  calculating  momentum  transfer  and  described 
in  Section  4.2.8.  Based  on  the  results  shown  in  Table  (D.l),  we 
found  that  the  major  source  of  energy  to  the  wind-generated  wave  is 
from  the  normal  pressure  p [column  heading  (1)  in  Table  (D.l)]. 

The  contributions  from  the  wave-induced  and  mean  turbulent  Reynolds 
stresses  [column  headings  (2),  (3)  and  (4)  in  Table  (D.l)]  are  ex- 
tremely small  (approximating  2.5%  of  the  total  energy  transferred  to 


the  wave  field).  The  surprising  result  is  that  the  terms  represented 
by  the  contribution  of  energy  to  the  wave  field  due  to  u'v'  and 
v'v'  are  negative.  This  implies  that  the  energy  is  extracted  from 
the  wave  field  and  transferred  to  the  wind.  This  result  can  not  be 
explained  at  the  present,  but  may  be  attributed  to  the  accuracy  of 
the  calculation  of  these  small  higher-order  terms.  As  indicated 
in  the  Appendix  A,  the  accuracy  of  calculation  of  such  terms  is  about 
30%.  But  according  to  the  reasonable  trend  of  the  results  of 
- pa  u'v'u,  - pa  v' 2v  and  pa  u'v'  9n/3x  v [column  headings  (2), 
(3)  and  (4)  in  Table  (D.l)]  with  increasing  free  stream  wind 
velocity  u„  and  fetch,  we  can  see  the  order  of  magnitude  of  these 
terms  is  reliable,  and  hence  the  contribution  of  these  terms  to  the 
energy  transfer  from  wind  to  waves  is  negligible.  Therefore,  we  can 
conclude  that  almost  all  the  energy  transfer  from  wind  to  waves  is 
due  to  the  normal  pressure  acting  on  water  surface. 

The  calculated  results  of  y„  , which  is  the  ratio  of  the  energy 

Ej 

transferred  into  waves  over  the  energy  transferred  into  water,  are 
tabulated  in  Table  (D.l)  for  different  wind  speeds  and  fetches.  It 
seems  that  is  independent  of  free  stream  wind  speed  and  fetch. 

The  calculated  average  value  of  y is  0.89  which  means  about  89% 
of  the  energy  transferred  into  water  goes  to  the  wave  field  directly 
and  the  remaining  11%  of  the  energy  goes  to  the  drift  current  field. 
This  analysis  is  based  on  the  assumption  of  the  drift  current  velo- 
city equals  3%  of  the  free  stream  wind  velocity  due  to  the  lack 
of  experimental  results  for  drift  current  velocity  in  our  channel. 

Based  upon  the  measured  pressure  pm(t)  and  the  waves  q(t), 
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the  spectral  density  of  energy  and  momentum  transfer  from  wind  to 


waves  through  the  normal  pressure  acting  on  wave  surface  can  be  cal- 
culated by  equations  (4.19)  and  (4.20),  respectively.  The  results 
are  shown  in  Figures  5.18  and  5.19.  Again  the  experimental  results 
are  similar  along  the  channel,  and  only  two  representative  sets  of 
data  at  a small  fetch  (6.51  m)  and  a longer  fetch  (15.66  m)  are 
shown.  Both  spectral  distributions  have  a shape  similar  to  that  of 
the  wave  spectrum  with  a peak  at  the  dominant  wave  frequency  and  a 
small  hump  at  the  second  harmonic.  Therefore,  most  of  the  energy 
and  momentum  transfer  occurred  around  the  dominant  frequency  of  the 
wave  spectrum.  The  magnitude  of  these  results  are  of  the  same  order 
as  Dobson's  (1971)  which  were  obtained  in  the  ocean. 

5. 5 Energy  Transfer  of  Nonlinear  Wave-Wave  Interactions 

5.5.1  Theoretical  Results  of  Nonlinear  Wave-Wave  Interactions 
Ihe  nonlinear  wave-wave  interaction  mechanism  proposed 
by  Hasselmann  (1962)  is  conservative  and  can  redistribute  energy 
among  various  frequencies  in  a given  wave  spectrum.  The  theoretical 
expression  of  the  rate  of  energy  transfer  due  to  this  nonlinear 
mechanism  is  a triple  integral  expression,  equation  (2.9).  As  dis- 
cussed previously,  Barnett's  parametric  equation  was  used  to  approxi- 
mate the  theoretical  results  for  nonlinear  wave-wave  interactions 
Sn^(f)  . Based  on  the  measured  wind  wave  spectra,  the  theoretical 
Snl(f)  calculated  from  equations  (2.17),  (2.18)  and  (2.19)  are  shown 
in  part  (a)  of  Figures  5.20  to  5.31.  This  series  of  figures  is 
presented  following  the  order  of  decreasing  wind  speed  and  decreasing 
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fetch.  The  spectral  distribution  of  Sn^(f)  consists  of  three  re- 
gions: a positive  and  relatively  low  frequency  region,  a negative 
and  intermediate  frequency  region,  and  a positive  and  relatively 
high  frequency  region.  This  implies  that  the  energy  transfers  are 
away  from  intermediate  frequency  region  to  the  lower  and  higher 
frequency  regions  in  a wave  spectrum. 

In  order  to  discuss  the  results  of  Sn^(f)  further,  the  results 
of  the  net  rate  energy  transfer  S(f)  are  presented.  The  spectral 
distribution  of  S(f)  is  similar  to  that  of  Sn^(f)  , namely,  posi- 
tive transfer  in  a relatively  low  frequency  region  and  negative 
transfer  in  the  intermediate  frequency  region  as  shown  in  part  (a) 
of  Figures  5.20  to  5.31.  The  results  for  S(f)  were  calculated  from 
equation  (2.5)  and  based  on  the  measured  wave  spectrum  E(f)  along 
the  fetch.  This  net  rate  energy  transfer  as  shown  in  equation  (2.5) 
is  the  sum  of  the  rate  of  energy  transfer  due  to  air- input  %n(f)» 
nonlinear  wave-wave  interactions  Sn^(f)  , and  white-cap  dissipation 
Scjg(f).  The  resulting  positive  transfer  for  S(f)  in  the  relatively 
low  frequency  region  and  negative  transfer  for  S(f)  in  the  inter- 
mediate frequency  region  indicate  a shifting  of  wave  energy  from 
intermediate  frequencies  to  the  low  frequency  region  along  the  direc- 
tion of  wave  propagation  as  shown  in  Figure  5.1.  The  spectral  dis- 
tribution of  S(f)  is  consistent  with  the  observed  wave  spectrum 
establishment  along  a fetch  for  constant  wind  or  at  a fixed  fetch 
with  increasing  wind.  Since  the  air-input  S^n(f)  and  white-cap 
dissipation  S(js(f)  are  purely  positive  and  negative  transfers, 
respectively,  in  all  frequency  regions  and  because  the  nonlinear 


to  S(f)  , the  nonlinear  wave-wave  interactions  mechanism  must  ba 
responsible  for  the  energy  shift  from  intermediate  frequencies  to 
low  frequencies.  Consequently,  this  stresses  the  importance  of 
Sni(f)  to  wave  growth  in  any  relevant  frequency  region. 

The  parametric  equation  results  of  S^Cf)  are  examined  in 
the  previously  mentioned  three  regimes.  in  the  relatively  low 
frequency  region,  the  peaks  for  both  positive  energy  transfer  of 
Sni(f)  and  S(f)  are  approximately  at  the  same  frequency,  except 
for  two  cases  of  short  fetch  (3.48  m to  6.51  m)  with  wind  speeds  of 
7.09  and  8.01  m/sec  as  shown  in  part  (a)  of  Figures  5.30  and  5.31, 
respective] y.  This  miss  match  in  frequency  implies  that  the  parame- 
tric equation  can  not  predict  Sn^(f)  well  at  the  short  fetch,  as 
Sj^Cf)  is  compared  with  the  measured  Snp(f)  . Although  Barnett 
(1966)  indicated  that  the  parametric  equation  is  developed  for  the 
fully  developed  sea  ; it  is  also  applicable  for  the  partially 
developed  sea  . Since  the  two  cases  are  in  the  very 
initial  stjage  of  the  wave  generating  zone,  the  fetches  are  too 
short  to  reach  the  applicable  range  of  the  parametric  equation.  For 
all  other  cases  the  parametric  equation  is  indeed  to  be  able  to  pre- 
dict Snj[(f)  sufficiently  well  in  the  relatively  low  frequency 
region . 

In  the  intermediate  frequency  region,  the  negative  peak  of 
Sn^(f)  and  the  positive  peak  of  S^n(f)  (or  ECf)  ) occur  at  the 
same  frequency  as  shown  in  Figures  5.20  to  5.31.  Since  the  net 
transfer  S(f)  is  the  algebric  sum  of  Sjn(f),  Snj(f)  and  S^^f) 


and  the  negative  transfer  of  SjgCf)  is  always  smaller  than  the 
positive  transfer  from  wind  to  waves  S^n(f)  as  shown  in  part  (b)  of 
Figures  5.20  to  5.31,  the  nonlinear  energy  transfer  Sni(f)  can  be 
attributed  largely  to  the  negative  transfer  of  S(f)  in  the  inter- 
mediate frequency  region.  Consequently,  the  energy  transferred  from 
wind  to  waves  in  a intermediate  frequency  region  is  always  transfer- 
red to  both  the  relatively  low  and  high  frequency  regions  throught 
the  mechanism  of  nonlinear  wave-wave  interactions.  In  general,  the 
parametric  equation  results  are  quite  adequate  and  acceptable  as 
pointed  out  by  Barnett  (1966)  and  Mitsuyasu  (1968b)  and  evidenced  by 
the  results  of  this  study. 

The  relatively  high  frequency  region  seems  to  begin  with  the 
frequency  at  which  the  net  energy  transfer  S(f)  equals  zero.  This 
implies  that  the  positive  energy  transfer  of  S-^n(f)  and  Snl(f)  is 
balanced  by  the  negative  transfer  of  Sds(f)  in  the  relatively  high 
frequency  region.  The  calculated  results  for  Snj^f),  based  on 
Barnett's  parametric  equation,  were  shown  in  part  (a)  of  Figures 
5.20  to  5.31.  The  calculated  values  of  Sn^(f)  are  too  large  to 
satisfy  the  conservative  requirement  of  the  nonlinear  wave-wave 
interaction  mechanism  (Mitsuyasu,  1968b).  To  avoid  overestimating 
Sni(f)  by  Barnett's  parametric  equation  in  the  relatively  high  fre- 
quency region,  one  can  use  the  conservative  nature  of  the  nonlinear 
wave-wave  interaction  mechanism.  A suggested  method  is  to  use  the 
parametric  equation  to  calculate  the  sum  of  the  energy  transfer  in 
the  relatively  low  frequency  and  intermediate  frequency  regions  where 
the  spectral  value  of  Sn^(f)  changes  from  negative  to  positive,  say 


at  a frequency  f . A net  negative  result  is  expected  due  to  the 
large  negative  transfer  in  the  intermediate  frequency  region.  In 
order  to  satisfy  the  conservative  law  for  Sn|(f)  over  all  fre- 
quencies, we  require  a positive  energy  transfer  equal  to  the  negative 
value  obtained  from  the  previous  calculation  to  be  spread  over  the 
frequency  range  between  fg0  and  the  cutoff  frequency  of  the  wave 
spectrum.  A smooth  decay  function,  such  as  exponential  decay,  could 
be  assumed  for  redistributing  the  positive  energy  transfer  in  the 
relatively  high  frequency  region.  With  modification  Barnett's  para- 
metric equation  would  be  the  most  applicable  method  for  predicting 
nonlinear  wave-wave  interaction  economically  and  efficiently. 

5.5.2  Experimental  Results  of  Nonlinear  Wave-Wave  Interactions 
Because  we  can  not  directly  measure  the  energy  transfer 
due  to  SjjjCf)  and  S^gCf)  separately,  the  experimental  result  of 
the  sum  of  Sn^(f)  and  Sds(f)  was  first  deduced  from  equation 
(4.21)  and  is  plotted  in  part  (a)  of  Figures  5.20  to  5.31.  As  shown 
in  the  same  graphs,  the  negative  value  of  (Sni  + S^g)  at  the  peak 
frequency  if  the  intermediate  frequency  region  is  larger  than  that 
of  the  parametric  equation  result  for  Sni(f)  as  expected,  except 
one  case  at  a fetch  of  3.46  m to  6.51  m and  wind  speed  of  8.88  m/sec 
(see  part  (a)  of  Figure  5.29).  In  order  to  evaluate  the  experi- 
mentally determined  Sn^(f),  the  white-capping  dissipation  Scjs(f) 
must  be  subtracted  from  (S^  + S^g) . Because  we  can  not  directly 
measure  S(jg(f),  a theoretical  model  of  white-capping  dissipation, 
Sds(f)  had  to  be  assumed  in  order  to  obtain  Sn^(f)  from  equation 


(4.22).  Therefore  S , (f)  was  estimated  from  the  measured  results 
for  S(f)  and  %n(f)  [from  equations  (2.5)  and  (4.19),  respectively] 
together  with  Hasselmann's  white-cap  dissipation  model  for  Sjg(f) 
[equation  (2.21)].  The  results  for  Sn^(f)  are  shown  in  part  (b)  of 
Figures  5.20  to  5.31  and  compared  there  with  the  theoretical  results 
for  Sn^(f)  computed  by  Barnett's  parametric  equation.  In  general, 
the  qualitative  comparisons  are  satisfactory  in  the  relatively  low 
and  intermediate  frequency  regions,  except  for  the  three  cases  where 
the  wind  wave  development  at  small  fetch  (3.46  m to  6.51  m)  is  at 
its  initial  stage  as  shown  in  Figures  5.27  to  5.31.  The  relatively 
high  frequency  components  of  the  wind-waves  at  small  fetch  and  the 
inadequate  representation  of  Snl(f)  with  the  Barnett's  parametric 
equations  as  discussed  previously  may  be  the  reasons  which  contribute 
to  the  disagreement  for  these  cases. 

For  the  cases  in  good  agreement,  both  the  theoretical  Sn i ( f ) 
and  experimental  Sn^(f)  generally  have  similar  positive  peak  fre- 
quencies, similar  spectral  distributions  and  almost  the  same  amount 
of  positive  energy  transfer  in  the  relatively  low  frequency  region; 
they  also  have  similar  negative  peak  frequencies  and  similar  spectral 
distributions  in  the  intermediate  frequency  region.  The  positive  and 
negative  peak  frequencies  of  Sn^(f)  are  almost  the  same  as  the 
positive  peak  frequency  of  S(f)  and  peak  frequency  of  E(f)  , res- 
pectively. 

As  to  a quantitative  comparison  of  the  results  of  ^^(f)  and 
Snl(f)  in  the  intermediate  frequency  region,  the  parametric  equation 
results  for  Sn^(f)  are  always  slightly  larger  (5%  to  40%)  than  the 
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experimental  Sn^(f)  in  the  neighborhood  of  the  negative  peak  fre- 
quency for  the  cases  with  the  fetch  larger  than  6.51  m as  shown  in 
part  (b)  of  Figures  5.20  to  5.28.  These  discrepancies  might  be 
caused  by  an  overestimation  by  Hasselmann's  white-cap  dissipation 
model  for  S^gCf)  as  used  in  equation  (4.22),  resulting  in  an 
underestimation  of  Sn^(f). 

The  comparisons  of  S^Cf)  and  Sni(f)  in  the  relatively  high 
frequency  region  are  unsatisfactory.  The  experimental  result  for 
Snl(f)  did  not  indicate  the  expected  positive  transfer,  but  rather  ap- 
proached zero.  It  is  recalled  that  the  "experimentally"  determined 
Sjii(f)  was  calculated  from  equation  (4.22)  and  was  based  on  the 
experimentally  determined  S ( f ) and  Sin(f)  and  Hasselmann's 
Sds(f).  Since  the  measured  S(f)  is  approximately  zero  in  this 
relatively  high  frequency  region  as  shown  in  the  figures,  Sn^(f) 

~ ®in(f)  “ Sj s ( f ) . Based  on  the  experimentally  determined  Sin(f), 
it  seems  that  underestimation  of  S£js(f)  in  that  frequency  region 
may  be  responsible  for  Sn^(f)  not  being  able  to  have  the  expected 
positive  values  of  transfer.  As  discussed  previously,  unless  an 
independent  experimental  determination  of  the  dissipation  function 
S^gCf)  can  made,  ^^(f)  can  not  be  satisfactorily  evaluated 
over  all  frequency  ranges. 

★ 

After  completion  of  our  work  and  as  we  were  writing  it  up, 

Ramamon j iarisoa  and  Coantic  (1976)  reported  that  laboratory  wind- 
generated water  waves  with  frequencies  higher  than  the  dominant  wave 
frequency  do  not  obey  the  classical  deep  water  dispersion  relation- 
ship; these  waves  propagate  at  the  phase  speed  of  the  dominant  wave 
in  the  spectrum.  Consequently,  the  calculated  value  of  S(f)  should 
be  altered  before  further  detailed  comparisons  at  high  frequencies 
can  be  made. 
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Based  on  the  comparisons  and  discussions  presented  in  Sections 
5.5.1  and  5.5.2,  we  may  conclude  that  Barnett's  parametric  equation 


I 


* 


A 


might  slightly  overestimate  the  negative  transfer  for  Sn^(f)  in 
the  intermediate  frequency  region  and  is  inadequate  in  the  relatively 
high  frequency  region  (it  grossly  overestimates  the  positive  trans- 
fer). Hasselmann's  white-capping  dissipation  model  might  slightly 
overestimate  the  negative  transfer  in  the  intermediate  frequency 
region  and  underestimate  the  negative  transfer  in  the  relatively  high 
frequency  region.  However,  Barnett's  parametric  equation  is  appli- 
cable for  calculating  the  energy  transfer  due  to  nonlinear  wave-wave 
interactions  in  the  relatively  low  frequency  and  intermediate  fre- 
quency regions  of  a wind  wave  spectrum  with  the  normalized  fetch  x 
greater  than  100. 

5 . 6 Energy  Dissipation  of  White  Capping 

5.6.1  Theoretical  Results  of  White  Capping 

The  theoretical  results  for  white  capping  energy  dissipa- 
tion function  (f)  , calculated  from  equation  (2.21),  are  also 

shown  in  part  (b)  of  Figures  (5.20)  to  (5.31).  The  constant  d in 
equation  (2.21)  is  associated  with  the  energy  dissipation  spectrum. 

The  values  of  d were  calculated  based  on  the  energy  balance  in  the 
high  frequency  equilibrium  range  [equation  (4.26)]  and  are  tabulated 
in  Table  (5.2)  for  each  case.  The  graphs  clearly  show  the  distribution 
of  the  spectral  density  S(jg(f)  appears  similar  to  that  of  the  wave 
spectrum  as  expected  from  equation  (2.21).  The  dissipation  function 
used  herein  has  a maximum  negative  value  at  the  peak  frequency  of  the 
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wave  spectrum  and  decreases  on  both  sides  of  the  peak  frequency, 
finally  approaching  zero  in  the  relatively  low  and  high  frequency 
regions.  This  theoretical  model  of  S^Cf)  was  first  checked  by 
Hasselmann  using  JONSWAP  data.  Fairly  good  agreement  was  found  and 
presented  by  Hasselmann  (1974) . Because  of  the  lack  of  direct  mea- 
surements of  the  air-input  term  S^n(f)  , an  assumed  S^n(f)  was 
used  in  equation  (4.23)  by  Hasselmann  for  calculating  the  experi- 
mental results  for  Scjs(f).  In  this  study,  S^n(f)  was  measured 
directly  to  give  a more  accurate  investigation  on  Hasselmann' s white- 
cap  dissipation  model. 

5.6.2  Experimental  Results  of  White  Capping 

As  discussed  previously,  the  experimental  results  always 
yield  the  sum  of  Scjs(f)  and  Sn^(f).  In  order  to  deduce  the 
"experimental"  results  for  S(jg(f)  as  indicated  in  equation  (4.23), 
an  assumed  Sn^(f)  form  has  to  be  taken.  The  experimental  results 
for  S^gff)  are  evaluated  based  on  the  measured  values  of  S(f) 
and  Sj„(f)  and  Barnett's  representation  of  the  term  of  S ,(f). 

The  resulting  S^s(f)  are  plotted  in  part  (b)  of  Figures  (5.20)  to 
(5.31)  and  compared  with  the  theoretical  (f).  The  comparison 

between  S^ff)  and  S^g(f)  are  in  general  satisfactory  in  the 
relatively  low  and  intermediate  frequency  regions  for  the  fetches 
larger  than  6.51  m,  except  for  the  cases  with  fetches  between  6.51  m 
and  9.48  m and  12.61  m and  15.66  m and  with  wind  speed  of  7.09  m/sec 
[i.e.  Figures  (5.28)  and  (5.22)],  where  the  negative  peak  frequencies 
of  Sjjg^O  deviate  from  those  of  S^g(f).  The  Aviation  is  probably 


caused  by  the  overestimation  of  S , (f)  based  on  the  parametric 

ni 

equations  in  the  negative  peak  frequency  region.  Consequently,  the 
resulting  value  of  S,  (f)  would  be  underestimated  based  on  equation 

u S 

(4.23).  The  agreements  between  S^g(f)  and  S^g(f)  are  poor  in  the 
relatively  high  frequency  region  for  all  cases.  The  inadequacy  of 

Barnett's  parametric  equation  in  the  relatively  high  frequency  range 

* 

may  be  responsible  for  the  poor  agreement  (Mltsuyasu,  1968b).  The 
agreements  are  poor  for  the  small  fetch,  3.46  m to  6.51  m,  for  all 
wind  speeds;  this  may  be  attributed  to  the  inadequacy  of  Barnett's 
parametric  presentation  for  Sn^(f)  at  smaH  fetch  where  the  rela- 
tively high  frequency  waves  exist.  However,  generally  the  good 
agreements  are  shown  in  the  relatively  low  and  intermediate  frequency 
regions  for  large  fetches.  From  the  practical  point  of  view,  Hassel- 
mann's  dissipative  model  due  to  white  capping  is  acceptable  and 
applicable  for  developing  a new  wave  prediction  scheme  for  the  nor- 
malized fetch  x larger  than  100. 

5 . 7 Mean  Momentum  Transfer  from  Wind  to  Waves 

The  mean  momentum  transferred  from  wind  to  water  in  a developing 
turbulent  boundary  layer  calculated  from  equations  (4.29),  (4.30)  and 
(4.32)  is  tabulated  in  Table  (C.l).  Based  on  these  results,  we  found 
that  about  65%  of  the  momentum  transferred  into  water  goes  to  the 
wave  field  directly,  while  the  remaining  35%  goes  to  currents.  The 
most  important  contributions  to  this  momentum  transfer  are  the  normal 
pressure  acting  along  the  wave  surface  via  p 9n/9x  and  the  wave- 
induced  turbulent  Reynolds  stresses  acting  via  p iP~u' 9n/9x  . The 

Si 


former  contributes  about  73%  of  the  total  momentum  transferred  into 
waves,  while  the  latter  gives  about  27%.  Although  the  largest  contri- 
bution of  momentum  to  wave  growth  comes  from  the  normal  pressure,  the 
contribution  from  turbulent  effects  coupled  with  the  water  waves  can 
not  neglected.  Therefore,  the  turbulent-wave  interaction  must  be 
taken  into  account  for  the  momentum  transfer  in  formulating  wind-wave 
generation  theory. 

Since  the  phase  relations  between  pressure  and  wave  are  stable 
and  organized  over  a frequency  range  around  the  dominant  frequency 
of  the  wave  spectrum  as  shown  in  Figure  5.15  (f),  there  is  an  effi- 
cient momentum  transfer  for  each  frequency  component.  On  the  other 
hand,  the  result  of  calculating  pa  u'u'  9n/8x  for  the  momentum 
contribution  to  the  wave  field  indicated  an  unstable  phase  relation- 
ship with  the  wave  height,  resulting  in  random  fluctuating  values 
varying  from  positive  to  negative  throughout  the  frequency  range. 

Since  the  quantity  u'u'  is  a wave-induced  quantity  and  closely 
■j  correlated  with  the  wave  spectrum,  it  yields  a net  positive  transfer 

to  the  wave  field. 

The  averaged  ratio  y of  momentum  transferred  into  waves  over 
the  total  momentum  transferred  into  water  is  about  0.65.  Since  the 
equation  used  (Appendix  C)  in  obtaining  the  present  results  is  de- 
rived in  a developing  turbulent  boundary  layer,  the  mean  pressure 
gradient  in  the  wave  tank  has  already  been  taken  into  account.  In 
fact,  the  mean  pressure  gradient  in  our  wave  tank  is  small.  It 

l 

makes  our  results  comparable  with  the  field  results.  As  shown  in 
Figure  5.15  (b) , the  value  of  y varied  from  0.53  to  0.75  along 
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the  fetch  and  seems  to  approach  a constant  value  of  approximately 
0.70  at  large  fetch.  It  may  be  inferred  then  that  approximately 
70%  of  momentum  transferred  across  the  interface  goes  directly  to 
the  wave  field  for  fully  developed  wind  waves. 

Dobson's  (1971)  expression  of  y for  the  fully  developed 
boundary  layer  is 


a 


I 


y = - [(p  i?>/c  ] 


/ [-  0a  (u  + u' ) (v  + v' ) ] 


2tt f ~ 3n  , , 

T p « ' 


p (uv  + U*  v'  )] 
a ‘ 


which  is  different  from  the  present  expression  derived  in  Appendix 
C,  viz.  , 


y = (P  |J  + 0a  u'u'  / (p  + pa  u'u'  - pa  u'v') 


Since  p is  equal  to  - - p and  p is  also  equal 

3x  g v 3t  3x 

to  - pa  uv  for  the  inviscid  larminar  flow  as  indicated  by  Phillips 
(1966),  therefore,  the  only  difference  in  the  above  two  expressions 


of  y is  the  appearance  of  pa  u'u'  in  the  present  y . The 

neglectlon  of  this  term  which  is  a result  of  the  interaction  between 
turbulent  and  wave-induced  fields  could  cause  the  overestimation  of 
Dobson's  y value  [average  value  of  y is  1.1). 

Instead  of  direct  measurement,  JONSWAP  (1973)  calculated  y 
value  based  on  the  momentum  balance  scheme  of  a wind  wave  spectrum. 
These  values  of  y varied  from  0.1  to  0.8  as  discussed  in 
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Section  1.1.  However,  the  Dobson  (1971)  and  JONSWAP  (1973)  field* 
data  gave  a y value  of  about  0.8  for  a normalized  fetch  in  the 

^ r 

range  of  x = 10°  . The  present  y results  varied  from  0.53  to 
0.75  and  seemed  to  approach  a constant  value  close  to  0.7  for  nor- 
malized fetches  larger  than  100.  We  believe  that  the  experimental 
results  of  Wu  (1968)  and  Barnett  (1971)  and  Phillips  estimate  (1966) 
are  inconsistent  with  the  present  data  and  the  field  data  of  Dobson's 
and  JONSWAP' s.  Because  our  experimental  conditions  were  better 
controlled  than  other  investigators'  and  our  calculated  y values 
are  based  on  an  equation  consistent  with  our  experimental  conditions, 
we  confidently  conclude  that  about  70%  of  the  total  mean  momentum 
transferred  across  the  air-water  interface  goes  directly  to  the  wave 
field  in  the  fetch  range  of  100  < x < 500,  with  a reasonable  assur- 
ance that  this  result  can  be  extrapolated  to  larger  x values. 


6.  CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  Conclusions 


Based  on  the  presentation,  comparison  and  discussion  of  the 
present  results  in  Chapter  5,  the  following  conclusions  can  be  drawn: 
(1)  The  present  experimental  data,  as  well  as  the  preexisting 
data,  indicate  that  the  nondimensional  peak  frequency  f , nondimen- 


siorial  total  energy  e , and  the  equilibrium  constant  a of  a wave 

0.33 


spectrum  depend  on  nondimensional  fetch  x , e.g.,  f ^ x 


t a x1-23  and  a “ x ° , where  n varies  from  0.22  to  0.65  according 


to  the  available  data.  Since  these  parameters  satisfy  the  requirement 


of  a self-similar  spectrum,  i.e.,  e ~ a f ^ , with  better  under- 

m 


standing  of  these  parameters  we  will  be  able  to  derive  an  appropriate 
wave  spectrum  in  terms  of  normalized  fetch  x for  the  practical 
applicat ions . 

(2)  That  part  of  the  turbulence  field  close  to  the  interface  and 
with  frequencies  in  the  range  of  the  wave  frequency  is  strongly  coup- 
led with  the  wave  motion.  Strong  coherence  and  stable  phase  relation- 
ships between  p , u and  v with  respect  to  the  wave  field  q are 
found.  The  phase  angle  of  p leading  the  wave  by  225°  is  identical 
to  that  found  in  Elliott's  (1972)  field  data.  This  phase  relationship 
between  wave  and  wave-induced  pressure  may  provide  useful  information 
for  theoretical  or  numerical  modeling  of  the  momentum  and  energy  trans- 
fer from  wind  to  waves  and,  hence,  for  prediction  of  the  growth  rate  of 
the  wave.  In  a frequency  range  greater  than  10  Hz,  the  effects  of 
waves  on  the  turbulent  structure  were  not  detected  at  all. 
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inertial  subrange  of  the  turbulent  energy  spectrum  satisfied  the 
-5/3 

Kolmogorov  f Law,  which  is  consistent  with  the  physical 

assumptions  of  the  dynamic  behavior  of  a turbulence  spectrum  in  the 
equilibrium  range.  In  the  high  frequency  range  (f  > 500  Hz)  the 
f ^-law  characterizing  the  usual  viscous  range  of  a turbulence  spec- 
trum was  not  achieved  in  these  experiments. 

(3)  The  wave-induced  pressure  appears  to  play  a dominant  role 
in  energy  transfer  from  wind  to  waves.  The  wave-induced  turbulent 
Reynolds  stresses  u'v'  and  v'v'  , as  well  as  the  mean  turbulent 
Reynolds  stress  u'v'  have  small  effects  on  the  energy  transfer. 

The  surprising  negative  mean  transfer  due  to  the  wave-induced  turbulent 
Reynolds  stresses  may  be  attributed  to  the  inability  of  our  measurement 
system  to  resolve  the  very  small  wave-induced  Reynolds  stress  effects. 
Since  the  wave-induced  Reynolds  stresses  are  associated  with  non-linear 
effects  and  the  present  set  data  was  taken  at  a fixed  distance  of 
0.508  cm  above  the  highest  wave,  experimental  data  obtained  in  the 
lmnediate  proximity  of  the  interface  by  employing  a wave  following 
device  may  shed  some  additional  light  on  these  processes. 

(4)  The  nonlinear  wave-wave  interaction  mechanism  proposed  by 
Hasselmann  (1962)  does  indeed  provide  the  needed  transfer  process  from 
the  intermediate  frequency  region  into  both  the  relatively  lower  and 
higher  frequency  regions.  The  transfer  from  the  intermediate  frequency 
to  the  lower  frequency  region  provides  a logical  explanation  for  the 
establishment  of  a wind-wave  spectrum,  namely,  the  energy  is  shifted 
toward  the  lower  frequency  region  during  generation.  On  the  other 
hand,  the  transfer  from  the  intermediate  frequency  to  higher  frequency 


:X; 


region  provides  an  added  source  of  energy  to  balance  the  wave  dissi- 
pation process.  The  results  indicating  a large  positive  transfer  in 
the  high  frequency  region  due  to  nonlinear  wave-wave  interaction  may 
be  attributed  to  the  inadequacy  of  Barnett's  parametric  presentation 
of  the  nonlinear  process  in  that  frequency  range.  However, 

Barnett's  parametric  equations  (2.17)  are  generally  correct  and 
applicable  for  non dimensional  fetch  x greater  than  100.  Therefore, 
they  are  useful,  then,  for  practical  wave  prediction. 

It  seems  to  be  reasonable  to  conclude  that  wind  wave  generation 
and  growth  follow  this  process:  Soon  after  a spectrum  develops,  the 

wave-wave  interaction  mechanism  takes  over  to  transfer  energy  from 
intermediate  frequencies  to  lower  and  higher  frequencies  with  the 
atmosphere  acting  as  a continuous  source  of  energy  at  the  intermediate 
frequencies. 

Positive  energy  transfer  in  the  low  frequency  region  is  respon- 
sible for  the  evolution  of  the  wave  spectrum  during  the  period  of 
wave  growth  and  positive  transfer  in  the  high  frequency  region  is 
balanced  by  the  white-capping  dissipation. 

(5)  Hasselmann's  (1974)  model  for  white-capping  dissipation 

seems  to  be  valid  for  x > 100.  However,  this  validity  is  associated 

with  the  simultaneous  application  of  Barnett's  parametric  equation 

representing  the  nonlinear  wave-wave  interaction.  An  independent 

^ validation  of  either  the  S . (f)  or  S,  (f)  must  await  further 

nX  as 

experimental  studies  in  which  Sn^(f)  and  S^a(f)  can  maaaure<* 
separately. 

(6)  The  mean  momentum  transfer,  in  the  direction  of  wave 
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propagation,  from  wind  to  waves  in  a developing  turbulent  boundary 
layer  is  governed  by  the  normal  pressure  and  the  wave-induced  turbu- 
lent Reynolds  stresses  acting  along  the  wave  surface.  The  former 
contributes  about  73%  of  the  total  mean  momentum  transferred  into 
waves,  whereas  the  latter  contributes  about  27%.  Accordingly,  turbu- 
lence-wave interaction  plays  an  important,  but  not  dominant,  role  in 
the  momentum  transfer  from  wind  to  waves.  About  70%  of  the  total 
mean  momentum  transferred  goes  from  the  wind  to  the  wave  field  directly 
(in  the  normalized  fetch  range:  100  < x < 500),  while  the  remaining 

30%  goes  to  the  current  field. 

6.2  Recommendations 

(1)  The  measurement  simultaneously  of  wave  height  and  the  fluc- 
tuating pressure  and  velocity  components  in  a fixed  or  in  a wave 
following  system  over  the  surface  of  a mechanically-generated  water  wave 
is  a logical  extension  from  the  present  experimental  program.  This  would 
produce  a greater  influence  by  organized  wave-induced  components  and 
hence  better  estimates  of  the  energy  and  momentum  transfer  from  wind  to 
waves  associated  with  wave-induced  or  wave-associated  components.  With 
these  results,  we  might  be  to  formulate  a quantitatively  correct  model 
for  the  momentum  and  energy  transfer  from  wind  to  waves.  In  addition, 
based  on  the  experimental  data  on  the  energy  transfer  from  wind  to 
waves,  nonlinear  effects  seem  to  be  important.  Consequently,  experi- 
mental data  in  the  close  proximity  of  the  interface  which  could  be 
obtained  with  a wave  following  system  without  damaging  the  hot-wire 
would  be  extremely  valuable. 
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(2)  Barnett's  parametric  equation  and  Hasselmann's  white-capping 
dissipation  model  might  need  to  be  modified  in  the  relatively  high 
frequency  region  in  order  to  match  the  experimental  results  and  balance 
the  energy  transfer  in  this  region.  For  the  complete  verification  of 
these  two  models  and  the  full  understanding  of  energy  transfer  of  wind- 
generated gravity-wave  spectrum,  an  equivalent  experimental  program 
to  the  present  one  should  be  carried  out  in  the  field. 
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APPENDIX  A 


UNCERTAINTY  ANALYSIS 


Uncertainty  in  measurements  always  exists  due  to  improper  per- 
formance, inaccuracy  of  instruments,  and/or  propagation  of  input 
uncertainties  through  the  data  reduction  procedure  to  the  results. 
Kline  and  McClintock  (1153)  proposed  a procedure  to  determine  the 
percentage  uncertainty  in  the  general  dependent  quantity  of  interest 
f(x^)  as  follows: 
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where  the  x^  are  the  measured  inputs  in  the  data  reduction  procedure 
and  6f  is  the  uncertainty  in  f . The  experimental  uncertainties 
computed  here  are  based  on  the  odds  of  20  : 1, 

The  basic  percentage  uncertainties  in  6x^/f  were  calculated 


from  the  calibration  results  and  are  listed  below: 

Calibration  of  crystal  pressure  transducer  ± 5% 
Calibration  of  hot  films  ± 5% 
Calibration  of  wave-height  gauge  ± 4% 
Mean  velocity  ± 4% 


The  uncertainty  arising  from  spectral  estimates  is  computed  by 
the  following  relationship  [Bendat  and  Piersol  (1971)]: 
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where  N is  the  number  of  spectra  averaged  and  e is  the  uncertainty 
of  the  computed  spectra.  Fifty  spectra  were  averaged  and  used  for  the 
results  presented  here  so  the  basic  uncertainties  of  these  spectra  are 
± 14.2%.  The  corresponding  degrees  of  freedom  [=  2N,  Bendat  and  Piersol 
(1971)]  and  confidence  band  [Blackman  and  Tukey  (1958)]  are  100  and  90% 
respectively. 

The  resulting  uncertainties  calculated  by  equation  (A.l)  are 
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± 10% 
± 15% 
± 15% 
± 15% 
± 14% 
± 14% 
± 14% 

± 14% 

± 30% 

± 30% 

± 30% 


The  uncertainty  here  does  not  invalidate  the  results  and  conclu- 
sions presented  in  Chapters  5 and  6. 
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APPENDIX  B 


ENERGY  TRANSFER  EQUATION  OF  A WIND  WAVE  SPECTRUM 

Hasselmann  (1968)  posed  a generalized  formulation  of  the  energy 
balance  equation  for  a wind  wave  spectrum;  this  equation  can  be  ex- 
pressed as 


DF(k<X>t)  .11  ^ !Z_+k  II sfk} 

Dt  9t  Xi  9Xl  i 3ki  ‘ 00 


(B.l) 


where 

• _ 9(o  (x,£) 

i 3k±  ’ 

It  is  assumed  that  the  bottom  depth  H is  slowly  varying,  i.e., 

(1/k)  (9H  / 3x^)  <<  1,  so  that  the  geometrical  refraction  law  equation 

(B.2)  is  applicable  (the  dot  denotes  differentiation  with  respect  to 

1/2 

t and  to  * (gk  tank  kH)  is  the  radian  frequency  of  a free  surface 

gravity  wave  traveling  in  water  of  depth  H) . S(ic)  is  the  net  rate  of 

energy  transfer  and  is  a function  of  the  wave  number  vector  ic.  Also 

f F(Ic,x,t)  die  * p g <r?  > where  p is  water  density,  q is  the  water 

“ 0)  (0 

surface  deviation  from  the  mean  still  water  level  and  the  bracket  de- 
notes an  ensemble  average. 

In  equation  (B.l),  DF  / Dt  represents  the  Lagrangian  rate  of 

change  of  the  spectrum  relative  to  a wave  group  ic  moving  along  the  path 

in  x - ic  phase  space  determined  by  the  Hamiltonian  equations  (B.2). 

• 

In  deep  water,  the  wave  refraction  k^  9F  / 9k^  is  zero.  The  source 


• 9m  (x,ic) 

i 3xa 


function  S(k)  represents  the  net  rate  of  energy  transfer  to  or  from 
the  spectrum  at  the  wave  number  & due  to  all  the  interaction  processes 
which  affect  the  component  1c.  The  general  form  of  the  source  function 


S(k)  - S.  + S.  + S,  + S.  + S.'  + S'  + S'  + S'  + Se  + S,  + 5,  + ... 
12341234567 


represents  the  constant  energy  transfer  to  the  wave  field 
through  atmospheric  turbulence  pressure  fluctuations  (Eckart, 
1953;  Phillips,  1957). 

represents  the  inviscid  mechanism  of  unstable  coupling  be- 
tween the  wave  field  and  the  mean  boundary  layer  flow  (Miles, 
1957). 

represents  the  non-linear  correction  to  Miles'  inviscid 
theory . 

represents  the  energy  transfer  due  to  wave-burbulence  inter- 
action (Hasselmann,  1967). 

represent  the  transfer  due  to  the  interactions  of  waves  with 
mean  currents  and  turbulence  in  the  ocean,  corresponding  to 
the  transfer  processes  of  similar  interactions  with 

the  atmosphere;  s|  and  are  normally  neglected  according 

to  Hasselmann  (1968). 

represents  the  energy  transfer  due  to  nonlinear  wave-wave 
Interactions  (Hasselmann,  1962). 
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represents  energy  dissipation  in  shallow  water  due  to  turbu- 
lent bottom  friction  (Hasselmann  and  Collins,  1968). 
represents  energy  loss  due  to  wave  breaking  (Hasselmann, 
1974). 

...  other  mechanisms  for  a systematic  and  complete  theory  of 
wind-wave  generation. 

Equation  (B.l)  can  be  simplified  for  a stationary  deep-water  wave 
spectrum  with  negligible  of  mean  ocean  current  and  can  be  expressed 
as 


Sin  + Snl  + Sds  (B.3) 

where  x is  the  group  velocity. 

Next,  we  let  Sin  represent  the  rate  of  energy  transfer  from  wind 
to  waves,  corresponding  to  the  sum  of  to  S^,  which  is  dominated 

by  the  normal  pressure  acting  on  the  wave  surface  as  indicated  in 
Appendix  D;  and  are  equal  to  S,.  and  S^,  respectively. 

Equation  (B.3)  is  equivalent  to  equation  (2.5)  in  Section  2.1. 
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APPENDIX  C 


i 


MOMENTUM  TRANSFER  FROM  WIND  TO  WAVES 
IN  A DEVELOPING  TURBULENT  BOUNDARY  LAYER 

The  developing  turbulent  boundary  layer  equation  describing  the  mo- 
mentum transfer  across  the  wind-wave  interface  was  derived  by  C.  T.  Hsu 
(private  communication)  and  is  summarized  in  this  section.  In  the 
derivation  of  the  governing  equation  of  the  air  flow,  the  decomposition 
of  an  instantaneous  quantity,  f(x,t)  for  the  flow  above  the  wavy 
boundary  is  assumed  to  be 

f(x,t)  - f (x)  + f (x , t ) + f'(x,t) 

where  f is  the  time-averaged  mean  quantity,  f is  a wave-induced 
perturbation  quantity  associated  with  the  wave  motion,  f'  is  a back- 
ground turbulent  quantity.  In  order  to  determine  the  wave-induced 
quantity,  the  statistical  techniques  of  time  average  and  phase  average 
are  employed.  The  time  average  of  f(x,t)  is 

f(x)  ■ lim  ^ J f ( x , t ) dt 
1 o 

whereas  a phase  average  is 

N 

<f(x,t)>  - lim  — f (x,  t + nr) 

N-k*>  N n 1 

where  t is  the  period  of  a wave  and  N is  the  number  of  waves  in  the 
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ensemble.  The  background  turbulence  is  assumed  to  be  random  and  its 
contribution  to  the  phase  average  of  a large  ensemble  is  zero.  Then 
the  phase  average  of  f is 


<f>  - f + f 


and  the  wave-induced  quantity  is 

i - <f>  - 7 


Some  useful  relations  that  follow  from  the  above  definitions  are 


f'  - 0 

f - 0 

<7>  - f 

7f7  . 

7 g - 7 g 

<f  g> 

f g’  - o 

A 

"(SO 

I'M 

V 

<f'>  - 0 
<f>  - f 

f<g>  <f  g>  - f<g> 

0 < f g'>  ■ 0 


where  f and  g are  two  arbitrary  quantities. 

C.l  Momentum  Transfer  in  Horizontal  X-Plrection 

The  air  flow  above  the  wave  surface  n is  assumed  to  be  a 
two-dimensional  turbulent  shear  flow.  The  momentum  equation  in  x- 
dlrectlon  and  continuity  equation  can  be  written  as 


(C.l) 


i!i+  u i H + viH 
3t  3x  ay 


1 3p  . /32u  , 3^u, 


|a  + |i  - 0 

3x  3y 


(C . 2) 


Applying  (C.2),  we  can  rewrite  (C.l)  as 


lit  + ini  + 2liUl  . _ _I  |£+vV2u 
3t  3x  3y  Pa 


(C.3) 


As  previously  discussed  we  can  decompose  the  velocity  components  and 
pressure  into  mean,  wave-induced  and  background  turbulent  quantities 
as  follows: 


u ■ u + u + u' 


v ■ v + v + v' 


(C.4) 


P - P + P + P' 


Substituting  (C.4)  into  equation  (C.l),  and  taking  the  phase  average  of 
the  result  yields 


1^-  ( TT+  G)  + [(u  + u)2  + <u,2>]  + |y  [(u  + u)  (v  + v)  + <u'v’>] 


~~  |r  (p  + p)  + v v2  (u  + u) 
pa  ox 


(C.5) 


where  ■ u'^  + u'^  and  <u'v'>  - u'v'  + u'v'.  The  terms 


>vr 


% 

i 


■* 

s 


associated  with  the  background  turbulence  disappear  in  the  phase 
average  except  for  the  quadratic  terms  u'2  and  u'v'.  For  the 
convenience  of  the  following  derivation,  the  pressure  and  velocity  com- 
ponents are  redefined  now  so 

u«u  + u v • v + v P“P‘00+P’+P 

where  p^  is  the  pressure  at  free  stream.  By  applying  the  new  defi- 
nition of  u,  v and  p,  we  can  rewrite  (C.5)  as 


It  + lb  (“2  + <u'2>>  + <uv  + <u'v'>) 


13'  — . 2 

“ " (P  * P„>  + v V u 


The  momentum  equation  at  free  stream  is 


(C.6) 


3u»  . 3 u_ 

3 1 u«  3x 


_i  ajt 

Pa  3x 


(C.7) 


wflere  Uoo  and  are  the  free  stream  velocity  and  pressure. 

Subtracting  equation  (C.6)  from  equation  (C.7)  and  applies  the  still 
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It  (u-  ■ u)  + (u»u  - u2)  + (u«  - u>  $r + It  [v(u»  - u>^ 


'p7It(p  + p)+v72u"It <u'2>  ■ It <u’v,>  = ° {c-8) 


Now,  define  the  displacement  thickness  ^ and  momentum  thickness  6, 


•i-J,  « - > ‘2  ■ (1  - 1H  * 


The  integrating  equation  (C.8)  with  respect  to  y from  n to  <=° 
and  applying  the  identity 


f iL.  dy  - [f  f dy]  + [f  ^-] 

Jn  3xx  y aXl  yJ  1 3x!Jy-n 


where  the  variable  x^  represents  either  x or  t allows  equation 
(C.8)  to  be  expressed  as 


h V +n  <“-2  s2> + ».  + k H 


[v  |u  . <u'v'>  + Ip|i.viu|3-  + <u.2>  ill] 

3y  p 3x  3x  3x  3x  y-n 


(C.9) 


In  arriving  at  equation  (C.9),  the  boundary  condition  at  the  interface. 


i.e.,  v*|y  + u"3x  *s  aPP^^e<^* 
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The  left  hand  side  of  equation  (C.9)  represents  the  total  rate  of 
horizontal  momentum  transfer  across  the  interface  from  the  air  field. 
The  right  hand  side  of  equation  (C.9)  represents  the  mechanism  of 
Interface  transfer.  Since  we  are  not  interested  in  the  instantaneous 
time  varying  field,  the  time-average  of  equation  (C.9)  was  taken  and 
the  resultant  mean  momentum  transferred  across  the  wind-wave  interface 
is 


M 


r 3u 

lp  3y  ” 


, - 3jl  3u  la  , 
+ P 37-^3 x + 


P u'2  1^-] 
a 3x  y=n 


(C.10) 


In  equation  (C.10),  the  terms  associated  with  the  wave  form  n and 
the  wave-induced  quantities  should  be  account  for  the  mean  momentum 
transferred  into  wave  field  directly,  i.e., 


M 

w 


P 


3u  3n 

3x  3x  pa 


2 la, 

3xJy“n 


The  remaining  terms  should  be  attributed  to  the  mean  momentum  trans- 
ferred into  current  as  expressed  by 


M 

c 


[y 


3u 

3y 


- p 


u'v'] 


y=n 


If  the  location  of  measurement  (n  + 6*)  is  outside  the  viscous 
sublayer  of  the  air  flowfleld,  then  the  viscous  effect  can  be 
neglected  , and 
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(C.ll) 


* 

•» 


M " [‘  p u'v’  + P S+  pa  u’2 


y-n+«, 


r ■'*  3n  f 2 ^ n 1 

M„  - Ip  al  + pa  u a?W«, 


(C.12) 


M - [-  p u'v’l 


y-n+6. 


(C.13) 


Accordingly,  the  ratio  of  the  mean  horizontal  momentum  transferred  into 
waves  over  the  total  mean  horizontal  momentum  transferred  across  the 
wind-wave  interface  is  given  by 


Y - M / M (C.14) 

w 

This  expression  seems  to  be  applicable  for  a general  air  flow 

(developing  and  fully-developed  turbulent  boundary  layers)  and  imposes 

no  limitation  on  the  wave  forms  (small  or  finite  amplitude).  The 

results  for  M , M and  Y are  tabulated  in  Table  (c.l)  for  various 
w c 

wind  and  wave  conditions.  Discussion  is  provided  in  Section  5.7. 

C.2  Momentum  Transfer  in  Vertical  Y-Pirectlon 

Similarly,  the  phase  average  of  the  momentum  equation  in  the 
y-direction  is 


3v 

3t 


+ u-F'+v|i  + |-  <u*v*  > + <v'v’> 


3v 

3x 


3y  3x 


3_ 

3y 


13  , r , ,32v  , 32vs 


(C.15) 
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where  u ■ u + u,  v - v + v and  p ■ p + p + p 

00 

Integrating  (c  .10)  with  respect  to  y from  n to  ® produces 

It  v dy]  + h (uv  + <u'v,> " v ^)dy]  + v” + ^ 

" ^ p + <v'v,>  - v fy  + v £ » ‘ <u'v,>  f^y-n  (C>16) 

On  the  left  hand  side  of  equation  (C.16),  the  first  term  represents 

the  momentum  transfer  in  the  control  volume  bounded  by  x = x , 

Y 

x ■ + ax,  y “ n and  y = <*>;  the  second  term  represents  the  net 

force  acting  on  each  edge  of  x = x^  and  x = Xy  + Ax;  and  the  third 

and  fourth  terms  represent  the  force  acting  on  the  edge  of  y = <=°. 

The  terms  on  the  right  hand  side  of  equation  (C.16)  represent  the  total 
force  supported  by  the  water  field.  The  time-average  of  equation 
(C.16)  is 


k (/_  (uv  + <u'v,>  - v dy]  + vl  + poo 


trP  + v’v’  - V £ - <u'v'>  I5-] 


ay 


3xJy-n 


(C.17) 


The  last  term  on  the  right  hand  side  of  above  equation  is  the  only 
term  transferring  vertical  momentum  to  the  wave.  Therefore,  the  verti- 
cal mean  momentum  transfer  from  wind  to  waves  is 


M 

yw 


I"  P. 


<u'vf> 


1ST 

8xJy»n+d* 


t-  P. 


In  1 

3x  y*n+5* 


(C. 18) 
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This  vertical  momentum  transferred  from  wind  to  waves  will  act 


together  with  the  horizontal  momentum  transferred  from  wind  to  waves 
and  form  the  circular  motion  of  the  wave. 


APPENDIX  D 


ENERGY  TRANSFER  FROM  WIND  TO  WAVES 
IN  A DEVELOPING  TURBULENT  BOUNDARY  LAYER 


D. 1 Energy  Transfer  In  Horizontal  X-Dlrection 

The  energy  transfer  across  the  wind-wave  interface  can  be  derived 
by  multiplying  equation  (C.8)  by  2u  and  is  expressed  as 


fr - »2> + k » - ;-> + h ‘v*  - “3> + h t<»:  - o v] 


j2  2, 
3y  " « 


■ 2 u k <"’u’>  - 2 u w <u’v,>  * ^ fe  ® + P - 2 * ‘f^2 


+ 2 * 17  <“  - 2 » <fi>2  + 2 v t (» £>  ■ o 


(D.l) 


where  u * u + u,  v ■ v + v and  p ° p + p + p , 

00 

We  define  Dx,  6^  , 62  and  6^  as 

°x  “|})2  + <^)2]  dy 

{i  ■ L (1  ■ S?  dy 


52  " K d - dy 

^ H Uqq  Uoo 


6 - / a-  (1  - ) dy 


n u 


oo  V1‘ 


Employing  the  boundary  condition  v -3n  / 3t  + u3n/3xt  and  integrating 


I 


equation  (D.l)  with  respect  to  y from  n to  °°  yields 


k V + 'ik*i+k<'Z  V ~2kf  u<u’u,>  dy 


- 2[u<u'u'>  + 2 f <u'u’>  dy  + 2[u<u'v’>] 

3xJy*n  J 9x  7 y=n 

Tt 

+ 2 J*"  <u’v,>  |S  dy  - T f (P  + Udy  “ pT  [(P  + P>  u 3^]y=n 


a n 


+ — f (p  + p)-r“dy- 
PaJn  3x  pa 


2 v [u  |H]  + 2 f-  [/“  v u dy] 

9y  y«n  9x  J n 9x  7 1 


+ 2 v [u  -JH  I1]  - 0 

9x  9x  y=n 


(D.2) 


Dividing  equation  (D.2)  by  a factor  of  2 and  rearranging  it  yields 


[v'7!iu+T_PT1u  + <u'u'>  u - v-^T^u  - <u'v’>u] 


9y  Pa  8x 


8x  3x 


, Dx  n°°  r i , 9U  . , I 9u.  , 1 »"  — 9u  , 

+ _ . [<u  „ > -U  + <u  v . dy  . _ (p  + p)  _ dy 


18  < 2 X J.  1 28  X J.  1 8 / 2 J.  . 

2 (u»  62>  +2  u»^  61  +7^  (U-  V 


k (vf^~  <U’U,>  - *-p>  « dyi 

1 a 


(D.  3) 


On  the  left  hand  side  of  equation  (D.3),  the  first  five  terms 
evaluated  at  y ■ n indicate  the  energy  transfer  from  wind  to  water 
due  to  horizontal  force;  the  sixth  term  denotes  the  large  scale  dissi- 
pation in  the  wind  field  and  is  considered  negligible;  the  seventh 
term  denotes  the  energy  sink  to  the  turbulent  field  due  to  the  strain 
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against  the  turbulent  stresses;  and  the  eighth  tens  represents  the 
work  done  by  the  pressure  force  p associated  with  volumetric  change 
in  x-direction.  The  net  horizontal  mean  energy  transfer  across  the 
wind  wave  interface  can  be  derived  by  taking  the  time-average  of 
equation  (D.3),  yielding 


r 3U  . 

3y  U^y-n 


+ --  [p  I11  u] 
pa  3x  y 


+ [<u’u'>  u] 

3x  y» 


y-n 


- [<u'v'>  u] 


y-n 


(D.  4) 


Applying  u - u + u [c.f.,  equation  (C.6)]  in  the  above  equation,  we 
obtain 


E 

x 


“‘Vn  + 


^ “ y 


+ ~ [P  T^] 

" 3xJy-r 


1_ 

Pa 


u 

y=n 


+ T“  [p  UI  + [<u'ur>  u]  + [<u’u'>  u] 
n — * 3x  y*ri  j l •iv  ■* 


Pa 


3x  y-n 


3x  y=n 


- V [~  Y1  u] 

3x  3x  y-n 


v [li  in  [lQ  111  qi 

l3x  3x  Jy-n  l3x  3x  Jy=n 


- [<u'v'>  u]  - [u'v'  ul 

y-n  y-n 


(D.5) 


From  equation  (D.5)  the  terms  associated  with  the  mean  velocity  at  the 
interface  should  account  for  the  energy  transfer  from  wind  to  current. 
At  the  interface,  the  mean  velocity  can  be  represented  by  the  velocity 
of  the  drift  current.  The  remaining  energy  should  be  transferred  into 
the  wave  field.  The  energy  transferred  into  wave  and  current  fields 
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are  expressed,  respectively,  as 


E 

xv 


[v  u + p U + <u'u’>  u 
3y  Pa  r 3x  3x 


3u  3n  - 

— \j  — *■  u — 

3x  3x 


3u  3n  - 

v T~  “ u - 

3x  3x 


u] 


y-n 


(D.6) 


E 

xc 


Pa 


JLpl* 
Pa  3x 


+ 


<u'u'> 


In 

3x 


- v 


3u  3 n 
3x  3x 


u'v’)  ul 

y-n 


(D.7) 


* 


D. 2 Energy  Transfer  in  Vertical  Y-Direction 

Similarly,  the  energy  transfer  equation  in  the  y-direction  can  be 
derived  by  multiplying  the  momentum  equation  (C.15)  by  2v  and  inte- 
grating it  with  respect  to  y from  n to  °°,  viz.. 


<v  Ty  » V,  - k '<*  + » v Vn  - I ■ ^ > "y-n  - "S  ^ " 


y-n 


+ [<u'v  > v]  + — ^ - r [<u  v >~  + ■cv’v’ >4^]  dy 

3x  y-n  Pa  J 3x  9yJ  y 

n 

- k1'  ♦ f)  - - **  v;  v2 «)  - 1 1 it;  ■ v2  w 


” 2 V®  “ L (<u'v,>  ■ V If ) V dy] 


3x 


3x 


(D.8) 


where  u - u + u,  v - v + v,  p-p  + p + p and  D « p I’'0  [(-^)2  + 

(>V,  , ‘ 7 " “ 

ty  1 dy 


y 
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Then,  the  vertical  mean  energy  transfer  across  the  wind-wave  interface 
can  be  derived  by  taking  the  time-average  of  equation  (D.8),  yielding 


Ey  " v ~ ( P + v “ pa  <v'v,>  v - V-tJ  g v + Pa<u’v,>  V]y,n 


(D.9) 


Substituting  v ■ v +v  « v . we  get 

n n n n 


Ey  " " P 0 " pa  <v'v,>  v - nggv  + pa  <-'▼'>  -g  v]^ 


(D. 10) 


Viewing  equation  (D.10),  we  find  that  all  the  vertical  mean 
energy  transfer  goes  directly  to  the  wave  field.  Finally,  we  can  sum 
up  the  horizontal  and  vertical  mean  energy  transfer  from  wind  to  waves 
as  follows: 


E » E + E 
w xw  yw 


[p  fy  “ “ pa  “'v'  u + (p  + pa  u'u*  * P g)  g u]y_n 


+ t11  v ’ P v - Pa  vST'  v + pa  u'v'  g v 


+ u'v'  g V - p g g v] 

3x  M 3x  3x  Jy-n 


(D. 11) 


After  neglecting  the  viscous  effect  and  only  considering  the  lower 
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order  terms,  we  can  rewrite  the  total  mean  energy  transfer  to  waves  as 


E„  - { - pa  u'v'  u - p v - pa  v'v'  v + pa  u'v'  v]y_^  (D.12) 


while  the  total  mean  energy  transfer  across  the  wind-wave  interface  is 


E-E  + E + E + E 

xw  yw  xc  yc 


- E + E 
w xc 


[ - pa  u'v'  u - p V - Pa  v'v'  V + P„  u'v'  -I*  v] 

“ ox  y^n 


+ Pa 


t(',v  + r-PT!+  w>#i-  via  in.  - 


3y  Pa  F 3 


a "a u’v')  if) 

3x  3x  y»n 


(D. 13) 


where  E - 0 
yc 

The  ratio  of  the  mean  energy  transferred  into  waves  over  the  total  mean 
energy  transferred  across  wind-wave  interface  is 


* 


*E  " Ew  ' E 


(D.14) 


The  calculated  results  for  each  term  in  equation  (D.12)  were  shown 
in  Table  (D.l)  from  our  experiments.  It  was  found  that  the  normal 
pressure  working  on  wave  surface  is  the  dominant  mechanism  for  energy 
transfer  from  wind  to  waves.  The  average  value  of  v is  about  0.9 
which  means  most  of  the  energy  transferred  across  the  wind-wave  inter- 
face goes  to  wave  field  directly. 
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APPENDIX  E 


DATA  REDUCTION  COMPUTER  PROGRAMS 

The  most  Important  computer  programs  used  in  data  reduction  are 
subroutines  for  the  cross  spectrum  and  calculation  of  nonlinear  wave- 
wave  interaction  by  Barnett's  parametric  equation.  Based  on  these 
two  programs,  we  can  calculate  the  power  spectra,  coherence,  phase 
relation,  energy  transfer  from  wind  to  waves,  nonlinear  wave-wave 
interactions,  and  white  capping  dissipation. 
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CCROSS  SPECTRUM  SUBROUTINE 

CNs512. NUMBER  OF  OAT*  POINT  FOR  EACH  SET  OF  DATA  IN  FFT 

CNNl=2*W-l  ,MsN/10,K=N/?M.L  = 2*N*3,TAPER  = COSIN£  TAPER 

CNN1 ,M,K,L, TaPERsALL  DEFINED  IN  MAIN  PROGRAM 

LX  1 1024) sSTORE  TWO  SETS  OF  DATA  FOR  FFT 

L**» TTsAUTO-SPECTRUM  OF  TWO  SETS  OF  DATA 

CRXY, TXVsCOSPECTRUM  AND  QUADRATURE  SPECTRUM  BETWEEN  TWO 

CSETS  OF  DATA 

SUBROUTINE  CROSSIN, NN1,M,*,L) 

DIMENSION  XU024)  , TAPER(Sl)  .XXC257)  . YY  (257)  ,RXY(257  ) , 

1 TXT  1257 1 

COMMON  RXY, TXT, XX, YY,X, TAPER 

NN=2*N 

SUMlsO. 

SUM?sO. 

00  40  1=1, NN 1,2 
SUM1=3UM1*X (I) 

40  SUM2sSUM2*X(I+l) 

SuMlsSUMl/FLOAT (N) 

SUM2sSUM?/FL0AT(N) 

00  45  l=l,NNl,2 
XllJsX (l)-SUMl 
45  X(IM)sX(I  + l)-SUM2 
DO  59  1=1, M 
11=2*1-1 
HlsNN-II 

X(I1)=X(IIJ*TAPERII) 

xcii*i)=x(ii*i)*rAPERm 

X(III)=X(III)*TaPER(I) 

58  XU1I*1)=X(III*1)*TAPER(I) 

CALL  FFT(X) 

XXll)=xX(l)*a.*X(l)**2 

TYCn=YYll)*4.*X(2)**2 

RXm)=RXY(l)*2.*x(l)*x(2) 

TXtmsO. 

00  122  1=2, K 
11*2*1-1 
1 1 I=L-2* I 
RX*X(II)*X(IIl) 

TUlKlItl  j-XUII  + 1) 

XX (I)sxx (I)*RX*RX+TX*TX 
RY*X(II*1)*X(1II*1) 

T Y*X (1 1 1 ) -X ( 1 1 ) 

YY(IJ=YY(IJ*RY*RY*TY*TY 

RXV(l)cRXY(l)*X(H)*X(lIl*l)*X(III)*X(II*l) 

122  TXYU)=TXYimx(lI)**2*Xim>l)**2-XlIlI)**2 
1-X ( 1 1 1*  1 ) **2 
RETURN 
END 

LIST  END  •*** 
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csubroutine  of  the  rate  of  energy  transfer  due  to 

CN0NL1NE AR  WAVE-WAVE  INTERACTIONS  CALCULATED  BY 
CBARNEET'S  PARAMETRIC  EQUATION 

CYY*INPUT  WAVE  HEIGHT  SPECTRUM,  DF sB ANDIV1 DTH  OF  SPTR 
SUBROUTINE  WNIICYY.DF) 

DIMENSION  YYC257) 

COMMON  A(257),B(257),SWN(257) 

SUMVSO. 

SUMYF  »0 • 

DO  1 1*1*257 
F*OF*FlOAT (1-1) 

SUMYeSUMY  + VYd) 

1 SuMYF*SUMYF*YY(I)*F*F 

e*sumy*df 

FO* (SUMYF  * OF /£ ) * * 0 .5 
Fl0*(F0*10.)**7*£*E/980.66**4 
B 1*7. 5*F 10* Cl .*2048./ (45. *3. 14159) ) 
A1*4.4*F10*E*F0*10.*3.*3. 14159/8. 

C*960 . 66*0F 
F53*0.53*F0 
F42*0.62*F0 
A ( 1 ) *0 • 

Bll)*0. 

SWW(1)*0. 

DO  2 1*2*257 
FcOF*FLOAT(1-1) 

F5*F-F53 
IF  (F5-0 . ) 3*  3*  4 

3 B ( I ) *0 . 

GO  TO  5 

4 B(I)*B1*F5*F5*F5/F*YV  II)*C 

5 FlaFO/F 

F2*l • -0 . 42*F 1 
F3*l.-Fl 

F4*-4.*F3*F3*0. 1*(F1**5) 

IF (F-F  42)6*6*7 

6 A C 1 J *0  . 

GO  TO  8 

7 ACI)*A1*F2*F2*F2*EXP(F6)*C 
6 SWW l I ) * A ( I ) -B ( I ) 

2 CONTINUE 
WRITCC6,9) 

4 FORMAT ( lHl ) 

WRITE(6*13)E*F0»A1*8! 

13  F0RMAT(iX,"Es",El0.4,"F0="*El0.4,"Al=H,El0.4, "Bis" 
1*610.4) 

RETURN 

END 

LIST  ENO  **** 
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ini 




B fc  T I il  iB—  ■ ■ >- 


Friction  velocity 


Notes:  (1),  (2),  (3)  and  (4)  corresponding  to  fetches  between  3.46  m and  6.51  m,  6.51  m and  9.48 
9.48  m and  12.61  m,  and  12.61  m and  15.66  m,  respectively. 
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energy  transfer  In  a developing  turbulent  boundary  layer 
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Figure  2.1  Comparison  of  the  rate  of  nonlinear  energy 

transfer  calculated  by  Hasselmann's  theoretical 
[S£i(f)J  and  Barnett's  parametric  [Sni(f)] 
equations  (after  Barnett,  1966) 
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Figure  3.5  Detailed  structure  of  crystal  pressure 
transducer  (after  Yu,  el  al.,  1973) 
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Figure  3.7 
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Calibration  of  the  magnitude  of 
crystal  pressure  transducer  output 
versus  frequency 
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Figure  3.8  Calibration  of  phase  lag  of  the  sensing 
system  of  crystal  pressure  transducer 
versus  the  variation  of  frequency 
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Figure  4.2  Mean  velocity  profiles  at  fetches  of  12.61  and  15.66  m; 

series  I,  II  and  III  corresponding  to  free  stream 
velocities  of  7.09,  8.01  and  8.88  m/sec,  measured  at 
fetch  of  12.61  m,  respectively 
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Figure  5.1  Power  spectra  of  wave  height  for  different  fetches  and  wind  speeds 


normalized  total  energy  (e)  of  the  wave  spectrum  versu: 
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Figure  5.3  Comparisons  of  present  and  previous  laboratory 
and  ocean  results  of  (a)  Phillips  equilibrium 
constant  (a) , (b)  normalized  total  energy  (e) , 

and  (c)  normalized  peak  frequency  (fm)  of  the 
wave  spectrum  versus  normalized  fetch  (x) 

(after  Hasselmann,  et  al. , 1973) 


Figure  5.3  For  legend  see  pp.  129 
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Normalized  power  spectra  of  u'  and  v'  at  fetch  of  15.66  m 
fcr  wind  speeds  of  (a)  7.26,  (b)  8.10  and  (c)  8.97  m/sec 


Normalised  power  spectra  of  u*  and  v'  for  different 
wind  speeds  at  fetches  of  (a)  6.51  and  (b)  15.66  m 


Normalized  power  spectra  of  u'  and  v*  at  different  fetches 
for  wind  speeds  of  (a)  7.09,  (b)  8.01  and  (c)  8.88  m/sec 
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Power  spectra  of  wave-induced  pressure  p and  velocity 
components  u and  v,  and  spectra  of  coherence  (R)  and 
phase  angle  (q)  between  p,  u,  v and  n at  fetch  of 
6.51  m and  wind  speed  of  8.62  m/sec 
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dominant  frequency,  and  the  ratio  y of  momentum  transferred  into  waves  over 
the  total  momentum  transferred  across  the  wind-wave  interface  versus  the 
normalized  fetch  x 
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Figure  5.16  Comparisons  of  the  present  (///)  and  previous 
results  of  phase  angle  between  p and  n 
(after  Dobson,  1971) 


-143- 


»*# 


9.0 


(a)  FREQUENCY  (Hi) 


(b)  FREQUENCY  (Hi) 


(c)  FREQUENCY  (Hi) 


Figure  5.19  Spectra  of  energy  and  momentum  transfer  from  wind  to  waves 
due  to  normal  pressure  acting  on  water  surface  at  fetch 
15.66  m for  wind  speeds  of  (a)  7.26,  (b)  8.10  and 

(c)  8.97  m/sec 


(b)  FREQUENCY  (Hi) 

5.20  (a)  Spectra  of  net  energy  transfer,  S(f) ; air-input,  Sin(f) ; nonlinear  wave-wave  interactions, 

Snl(f)  [Samett];  and  experimental  results  of  S^Cf)  + $ds(f).  (b)  Comparisons  between  theoretical 

lSni(f),  Sds(f)]  and  experimental  [Snl(f),  Sds(f)]  results  of  nonlinear  wave-wave  interactions 
and  white-capping  dissipation  at  fetches  between  12.61  m and  15.66  m and  wind  speed  of  8.88  m/sec 


Figure  5.21  (a)  Spectra  of  net  energy  transfer,  S(f);  air-input,  Sin(f) ; nonlinear  wave-wave  interactions, 

Sni(f)  [Barnett];  and  experimental  results  of_  Sni(f)  + Sds(f).  (b)  Comparisons  between  theoretical 
[Sni  (f)  , S(js(f)]  and  experimental  [§nl(f) , Sds(f)]  results  of  nonlinear  wave-wave  interactions 
and  W'hite— capping  dissipation  at  fetches  between  12.61  m and  15,66  m and  wind  speed  of  8.01  m/sec 
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(a)  Spectra  of  net  energy  transfer,  S(f);  air-input,  S^n(f);  nonlinear  wave-wave  Interactions, 

Snl (f > [Barnett] ; and  experimental  results  of^  §ni(f)  + Sds(f) . (b)  Comparisons  between  theoretical 

[SnlCO,  Sdg(f)]  and  experimental  [^^(f),  §cjg(f)]  results  of  nonlinear  wave-wave  interactions 
and  white-capping  dissipation  at  fetches  between  12.61  m and  15.66  m and  wind  speed  of  7.09  m/sec 
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Spectra  of  net  energy  transfer,  S(f);  air-input,  S<n(f);  nonlinear  wave-wave  interactions, 
f)  [Barnett]  ;and  experimental  results  of  &nl(f)  + §ds(f)*  (b)  Comparisons  between  theoretical 

(f) » Sds(f)]  and  experimental  tSni(f),  Sds(f)]  results  of  nonlinear  wave-wave  interactions 
white— capping  dissipation  at  fetches  between  9.48  m and  12.61  m and  wind  speed  of  8.88  m/sec 
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Spectra  of  net  energy  transfer,  S(f);  air-inpat,  Sin(f);  nonlinear  wave-wave  interactions, 
f)  [Barnett];  and  experimental  results  of  SniCf)  + §ds(f)-  (b)  Comparisons  between  theoretical 

(f) » S<js(f)]  and  experimental  [§nl(f),  J>ds(f)]  results  of  nonlinear  wave-wave  interactions 
white-capping  dissipation  at  fetches  between  9.48  m and  12.61  m and  wind  speed  of  8.01  m/sec 
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Spectra  of  net  energy  transfer,  S(f);  air-input,  Sin(f);  nonlinear  wave-wave  interactions, 
f)  [Barnett];  and  experimental  results  of  Snl(f)  + Sds(f).  (b)  Comparisons  between  theoretical 
(f),  Sds(f)]  and  experimental  [§nl(f),  Sds(f)]  results  of  nonlinear  wave-wave  interactions 
white-capping  dissipation  at  fetches  between  3.46  m and  6.51  m and  wind  speed  of  8.01  m/sec 
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Spectra  of  net  energy  transfer,  S(f);  air-input,  Sj^Cf) ; nonlinear  wave-wave  interactions, 
f)  [Barnett];  and  experimental  results  of^  Sni(f)  + Sds(f).  (b)  Comparisons  between  theoretical 
(f)  , SdsCOJ  and  experimental  fSni(f),  S<js(f)]  results  of  nonlinear  wave-wave  interactions 
white-capping  dissipation  at  fetches  between  3.46  m and  6.51  m and  wind  speed  of  7.09  m/sec 


